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ABSTRACT 


The research area of this thesis is an operational issue of the Venezuelan Navy. 
The specific area of research is the development of a Maritime Traffic Control Model that 
guarantees the efficient accomplishment of surveillance and protection of the territorial 
sea. The principal task is reducing the number of potential illegal shipments, for example 
drugs, toxic waste, or other outlaw activities in the Venezuelan sea. 

Currently, several operational activities are executed as a response against illegal 
shipments. However, these operational activities require many resources and a 
considerable amount of time. These operational tasks depend primarily on intelligence 
efforts that represent a high financial cost and additional risky actions. For these reasons, 
the successful execution of the maritime control mission requires more dynamic and 
efficient approaches to maximize operational benefits. 

One solution for this problem is the development of a stochastic decision making 
model, to analyze and set up inspection or interdiction operations in those areas whose 
geographic features represent closed transit areas for navigation. This decision model 
should be an important aid for the execution of maritime traffic control operations in 


closed maritime areas in the Caribbean Sea under the sovereignty of Venezuela. 
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EXECUTIVE SUMMARY 

The research area of this thesis is an operational issue of the Venezuelan Navy. 
The specific area of research is to develop a Maritime Traffic Control Model that 
guarantees the efficient accomplishment of surveillance and protection of the Territorial 
Sea. The principal task is reducing the number of potential illegal shipments, for example 
drugs, toxic waste, or other outlaw activities in the Venezuelan sea. 

Currently, several operational activities are executed as a response against this 
threat. However, these operational activities require many resources and a considerable 
amount of time with negative tradeoffs. These operational tasks depend primarily on 
Intelligence efforts that represent a high financial cost. For these reasons the successful 
execution of the Mantime Control Mission requires more dynamic and efficient 
approaches to maximize operational benefits. 

One solution for this problem is the development of a stochastic decision making 
model, to analyze and set up inspection or interdiction operations in those areas whose 
geographic features represent closed transit areas for navigation. This Decision model 
should be an important aid for the execution of Maritime Traffic Control Operations in 
closed mantime areas in the Canbbean Sea under the sovereignty of Venezuela. 

The main mission of the Venezuelan Navy, during peace, is to guarantee the 
observance of the national laws and international rules inside its terntorial waters. One 
way to accomplish this mission is the permanent surveillance and traffic control for vessels 
in its waters. Venezuela is the neighbor of several countnes in the Canbbean Sea. By its 
geographical position, Venezuela is one of the principal mantime routes and sometimes 
serves aS an intermediate ‘base’ for the entire mantime traffic in South America. Also, 
because of its tropical and calm weather, Venezuela has a heavy tounst manitime traffic, 
mainly from other Caribbean colonies or islands. All these factors create a maritime 
environment that is often difficult to control and to enforce the observance of laws. 
Currently, illegal activities such as drug dealing, transport of chemical waste and other 
illegal activities are some of the new persistent problems. Their presence requires more 


forces and a different approach from that of conventional warfare to defeat them. 


This view makes such operations other than warfare missions a daunting duty. To 
face this problem, many questions arise: 1) What is the approximate number of ships that 
come in or come out from Venezuelan harbors? 2) What is the area with the biggest 
maritime traffic? 3) Given the current number of available resources, how we should 
allocate them to execute the desired maritime control in an efficient way? 4) Using the 
geographical characteristics of certain coast areas (Natural Navigation Channels), is it 
possibly to carry out control policies and inspect the maritime traffic, according to the 
current “intelligence” functions? 

This thesis develops the required decision model to assist in the process of 
allocation of inspection assets to areas. 

This study formulates a decision model with stochastic events. Its application 
exploits the geographical aspect of the Venezuela Coast and the current maritime scenario 
with potential threats like drug dealers and other illegal commercial activities. A model 
provides an analytical procedure to determine the best allocation of assets for the 
execution of a maritime control mission in closed areas. Asset allocation is based on the 
fraction of ships that are carrying illegal cargo in each area j during a time period t. This 
model uses forecast information or “intelligence” to deploy such assets according to a 
Bayesian estimation of the illegal fraction. Few quantitative data related to this issue are 
available, which makes difficult the use of classical statistical methods. 

This thesis should form the introduction document for future studies in this type of 
operational requirement, using a probabilistic approach and a decision model to measure 
the results of the use of policies. 

The model estimates the proportion or fraction of ships carrying illegal cargo 
which heads to harbor through natural straits or channels. Using the available 
‘{ntelligence” as input to generate an initial estimate of this fraction, the model formulates 
a decision making process as a guide to the planner for asset allocation. The initial 
allocation decision is based upon a weighting using the initial estimated mean of the 
fraction of illegal ships with illegal cargo, u;. Those areas to be inspected with higher 
fractions will receive more assets. The model demonstrates that if the decision maker 


decides to get a precise estimation of u, by carrying out several inspection with a fixed 


number of assets in each area, then he/she can maximize the number of illegal ships found 
with little change in the asset allocation. On the contrary, when the planner decides to 
make inspections only with the prior information from intelligence, he/she may reallocate 
assets more frequently, and more inspections are needed to get an estimated |1; value with 
a desired precision. Even though this thesis keeps separate the logistics considerations 
related to the assets allocation, the fact remains that reducing the uncertainty of the u, 
value should be helpful in minimizing the complexity of the allocation problem. 

This thesis involves the intelligence as a quantitative input into the decision making 
process. The Bayesian approach makes possible the use of “Subjective” information into 
the analytical model. The model is used to demonstrate that ‘intelligence’ can be 


translated as an important estimator in this stochastic analysis. 


Bolt 





I. INTRODUCTION 
A. BACKGROUND AND THEORETICAL FRAMEWORK 


This purpose of this study 1s to improve current Venezuelan Navy maritime control 
policies. In the accomplishment of its mission, the Venezuelan Navy is updating and 
creating new procedures and tactics to minimize the increasing number of illegal shipping 
activities in the Venezuelan territorial sea and interior waters. Also, financial and war 
material constraints require that more efficient methods are needed to execute this type of 


task and to ensure the achievement of peace and observance of the National Laws. 


B. CURRENT SITUATION 


The current Venezuela Maritime Control Doctrine is designed mainly to protect 
and to place under surveillance the Maritime traffic in a potential war situation. Though 
some procedures may be applied during peace situations, the doctrine’s war application 
uses navigational restrictions. The current day-to-day maritime threats occur in open and 
unrestricted navigation areas. This means expeditious analysis and more efficient plans of 
actions are required. These plans demand a large amount of information related to the 
maritime situation in certain areas or regions. Unfortunately, this information or 
‘{ntelligence” should be updated very quickly, given the interaction of lawless actions and 
the dynamic nature of the Caribbean area. Therefore operational procedures need to be 
highly flexible and quickly set, due to the temporal value of the “Intelligence.” The value 
of intelligence decreases rapidly with time. 

The contemporary threats require an immediate response. For example: an 
increasing and big menace is drug trafficking. Drug traffickers use every kind of craft to 
transport drugs. They use any route at anytime to guarantee the successful arrival of illegal 
goods. Their dangerous procedures and tactics require large numbers of maritime control 


resources. 


A constantly changing situation requires the ability to make swift decisions. Every 
day, dozens of ships navigate in Venezuelan waters. The Venezuelan Maritime area 
includes a large coast with many accessible sites. Several Sovereign State Islands lie close 
to Venezuela in their own territorial seas. This fact makes surveillance more difficult due 
to legal restrictions. A potential tactic of using natural channels (straits, rivers, river 
mouths) might simplify the development of ship inspection plans. These natural channels 
usually create “tows” of ships coming into or out of Venezuelan harbors. Also, Venezuela 
owns several small islands, which form a set of "narrow corridors." The use of focal patrol 
groups in these areas may allow a more efficient way to concentrate resources with 


subsequent favorable results. This situation is shown in Figure 1. 
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Figure 1. A ship inspection plan in a 
Natural channel 


This thesis develops the required decision model to assist in allocating patrol 


groups to areas to carry out ship inspections. 


C. THESIS GOALS AND OUTLINE 


This thesis formulates a decision model with stochastic events. The model exploits 
the geographical features of the Venezuelan coast and the current maritime scenario with 
its potential threats of drug dealers and other illegal commercial activities. This study 
responds to an operational requirement to support the Venezuelan Navy mission, which 
includes the enforcement of the national laws and the acquisition of the available means to 
provide enough capacity for this type of task. 

This thesis is intended to form the introduction document for future studies of this 
type of operational requirement, using a probabilistic approach to decision models to 
assess the results of potential policies. 

The model will have as input the number of ships traversing different natural 
channels during a time period, an assessment of the fraction of these ships carrying illegal 
cargo, the number of ships a patrol can inspect during a time period and an objective 
function. The output of the model will be a policy to allocate the finite number of patrols 
to the different natural channels. 

Unfortunately no available data exist as to the current number of ships that come 
into or out of Venezuelan harbors. Many outsider recreational and fishing ships travel in 
these waters. They usually do not report their presence to the local authorities. However, 
tankers or merchant ships are regulated in a better way. However, their number is much 
smaller than the number of small ships described above. 

In spite of the lack of data, a simulation can be used to assess measures of 
effectiveness for different policies. The decision model will specify the resource allocation, 


given its objective and the corresponding measures of effectiveness. 


D. THE MARITIME CONTROL TRAFFIC PROBLEM IN’ THE 
CARIBBEAN SCENARIO 


Currently the Venezuelan Navy executes maritime traffic control following 
traditional tactical procedures and deploys its forces in different areas along its coast line. 
The number of deployed ships corresponds to the availability of ships and is fixed. 
Sometimes, the planner is required to move ships from one area to another based on 
operational requirements, and the area losing the ships must sometimes wait for a long 
period of time to receive additional resources. Often during this time, many opportunities 
to intercept or deter illegal actions are missed, with substantial waste of time and 
resources. Further, constant operational activities compel tight execution schedules. The 
lack of a specified procedure to execute surveillance operations can divert resources 
toward activities more related to traditional operational tasks. 

Nowadays, Venezuela and other nations have signed several agreements to control 
and coordinate efforts against illegal activities such as contraband and drug trafficking. 
The importance and exigency of this current issue implies the implementation of a tactical 
tool to support efforts against illegal activities, to monitor the observed results and to 


make the necessary decision in “real time”. 


Il. FORMULATION OF THE DECISION MAKING MODEL 


A. MODEL ASSUMPTIONS 


te Intelligence 


A military operation should not be executed without prior knowledge of the 
probable scenario. To ensure success an efficient employment of resource information or 
‘{ntelligence” is required to decide the best action, given the objective of the mission. In 
this study, it is assumed some intelligence is available to forecast the presence of illegal 


shipments in a given area or scenario. 


2. Inspection Clearance 


Any ship in or entering Venezuela waters can be examined. This assumption 1s 


stated to ensure the randomness of the selection of a ship for inspection. 


3. Inspection Time and illegal shipment detection 


To make the model formulation more tractable in this thesis, it is assumed that the 
time to inspect a ship is known and the detection of an illegal shipment onboard an 
inspected ship is certain. This is the most important assumption of this model because of 


its impact on the randomness of the inspection process and complexity of the analysis. 


4, Number of required resources to inspect a ship 


It is assumed that the minimum number of required resources to inspect a ship 1s 


two vessels. This is based upon prior experiences regarding safety and backup. 


5. Asset’s inspection capacity 


A certain fixed asset’s inspection capacity is considered for this analysis. Given 
that the inspection time is known and fixed, it is possible to calculate how many ships an 


asset can inspect during a given period of length / 


6. Time to reallocate resources 


One important issue is the time to reallocate resources, after an inspection mission 
is executed and the decision maker decides to move resources from area j to area k. This 
time is considered to be negligible; it is assumed there is enough time between the time of 


the decision to reallocate and the start of the new inspection mission. 


qe Resources to be reallocated 


Regardless of its initial position, any unit or ship can be reallocated to a new region 
or area. No logistic or operational constraints affect the desired resource’s performance. 
Further, only an even number of ships will be repositioned due to the assumption that two 


ships are needed to make an inspection. 


8. Ships to be inspected 


For this investigation, it is assumed the probability of carrying an illegal cargo 1s 
the same for any ship randomly chosen in area j during a time period t. However, it is 
both possible and likely that the available ‘intelligence” suggests that some ships are more 
likely to be carrying illegal cargo than others. This consideration can be incorporated in 
further ennchments of this model. For example, a method for ranking the incoming ships 
and estimating the probability that a particular ship is transporting an illegal cargo could be 
created. Such a procedure would assist in determining which ship will be inspected and to 


maximize the expected number of illegal shipments detected. 


oS: Ships Arrival 


For this study, the ships arrival process is assumed as follows: each ship heads to 
harbor independently the others. The number of ships that has arrived to harbor during a 
time period (t, t+1) is independent of the number of incoming ships in a time period (t-1,t). 
Further, the rate of incoming ships does not depend on the time of year or specific 
activities during the year; it is considered constant. The arrival process is characterized by 


a flow rate, 4;; the expected number of incoming ships during a period length /. 


B. MODEL STRUCTURE 


This model represents a periodic decision making process. For example, each 
period may consist of hours, days, a week, a month, etc. At the start of each period, the 
planner decides how many patrol assets he/she will employ and how they should be 
allocated to different areas. The allocation determines the maximum number of ships that 
can be inspected in each area. At the end of each period, the decision maker analyzes the 
results and decides to move assets from one area to another or to keep the same asset 
distribution. Figure 2 is used to clarify the time notation used in this study. Period t refers 


to the interval (t-1, t]. 


Period Length 





Figure 2. Time representation for a Maritime Control 
M ission 


Inspection periods are of equal length 4 and periods are numbered sequentially 


starting at one. The index 0 refers to the starting information that is known before the first 


inspection starts. 


analysis: 


t 
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A; (t) 


Cap, 


dj(t) 


k(t) 
i j(t) 


Notation 


The following notation is used to describe the elements or factors involved in the 


index of time (see Figure 2). 
period length (hours). (see Figure 2). 

time to inspect a ship.(hours). 

number of areas to be patrolled, indexed j = 1,2,3....m. 

Flow rate of outside ships coming to harbor through area j during a 
time period of length / 

Fraction of ships entering area j that are “bad” (have an illegal 
shipment on board). This fraction is unknown and so Is treated as a 
random variable. 

Total number of assets that are available for assignment. 

number of assets assigned to area j in period t. An asset consists of 
a pair of ships. This is a decision variable. Also, we must have that 
5 Ait) <n for every t. 

= 

The number of ships that an asset can inspect per time period in 
area J. 

number of ships to be inspected in area j in period t. dj(t)max = 
Cap,A\(t). 

number of “bad” ships found in area } in period t. 


Updated mean of the proportion of ships in area j that are carrying 


illegal cargo after all inspections are carried out in period t. 


Go (t) Updated variance of the proportion of ships in area } that are 


carrying illegal cargo after all inspections are carried out in period t. 


2: The Objective Functions 


In a real scenario it 1s difficult to inspect every ship coming to port. First, to 
inspect a ship at sea is a very hard task with many risks, and is usually highly time 
consuming. Second, to inspect all ships would demand a huge number of resources and 
many hours of operations. These considerations are just some of the constraints present in 
the execution of such missions. 

Therefore, the problem is how to choose the number of assets, Aj(t), to be 
employed in area j during an inspection period t, to achieve one or more of the following 
objectives: 

A) Reduce the uncertainty on Xj. 

B) Maximize the expected number of “bad” ships inspected in area j, given the 
limit on inspection resources. 

C) Reduce the expected value of X; through deterrance. 

In analyzing these goals, there is an order of priority. The goal A) must come first 
because the proportion of illegal ships in each area must be known as accurately as 
possible before goals B) or C) can be accomplished. This priority order is followed. It 1s 
assumed that there 1s no change in the habitual ‘tactics’ or shipments procedures executed 
by “bad” ships during successive inspection periods. Once goal A 1s achieved, the decision 
maker can decide how many assets should be assigned in each area j to maximize the 
expected number of “bad” ships intercepted (Objective B). Later, the correct asset 
allocation will induce the reduction of the expected value of X;, due to the permanent law 


enforcement through deterrance.(Objective C). 


3. The Decision Algorithm 


In the planning of resources allocation, there 1s an “intelligence” input that predicts 
the percentage or fraction of ships that carry illegal cargo into each area j. This 


information is modeled using a Beta random variable, whose general pdf is : 


TEH+BO) ae 
T(a(t))*T (Bt) 


The steps of the procedure are as follows: 


f(p; (a(t), B(t)) ) = Pp) a ion ose 


(1) The fraction of ships entering area j in a time period that contains illegal cargo 
is assumed to have a Beta distribution with parameters o;(0) and 40). The intelligence 
information 1s assumed to provide: 

1 (0) = initial estimate of the mean fraction of ships carrying illegal shipments through 
area J, 

o *,(0) = initial estimate of the variance of the fraction of ships carrying illegal shipments 
through area Jj. 


The mean of a beta distribution with parameters o and B 1s 





a 
= , 2m 
Oat (2.1) 
and its variance 1s 
7 (ho a (2.2) 
(a+ B)°*(a+ B+1) 
Thus, the initial parameter values a,(0) and £(0) can be obtained as follows: 
One pe s(O0)* CO Eo) = 14(0)* 07 (0) (2.3) 
Oo (0) 
BO) 2 SE en (2.4) 


14(0) 
These are the parameters of the prior beta distribution for area j, given the forecast or 


‘intelligence’. | Equations (2.3) and (2.4) are used frequently in the remainder of the 


thesis. When used for some period t, the “0”’s in these equations are replaced by “t’’’s. 


(ii) Having determined @(0) and £(0), the planner chooses the number of assets 


to be deployed for the first inspection period in each area j, Aj(1). The number of assets 


deployed is determined by the following relationship: 


Aipsi(9) 
Cap) 
A; (1) = n} ——*—— | 72] = 1,2,...,m. ZS 
a1 = Cap 


Aj (1) is the closest integer if (2.5) does not result in an integer. 

(iu) After completion of the first inspection mission, the planner obtains the 
numbers of “bad” ships found in the first period, {kj(1)}. It 1s assumed that each kj(1) is a 
sample from a binomial distribution with parameters dj(1), number of ships inspected in 
area j, and probability of success, 44(0). This “new” information is employed to update 
the beta distribution parameters for every area. The “posterior” distribution is also beta 
with parameters: 

a(1) = @(0) + ki), (2.6.a) 
BQ) = B(0) + dQ) — ki). (2.6.b) 
From these, the updated values of the mean and variance of the proportion of “bad” ships 
in region j, 44(1) and o7,(1), are determined using Equations (2.1) and (2.2). 

(iv) A test is made to see whether to end the algorithm or to return to step (11). The 

decision criterion used will depend on the objective function and is discussed in the next 


section. 


Gs ASSETS ALLOCATION AND DECISION CRITERION 


The decision making cycle is repeated until the decision maker’s established 
objective has been achieved. Depending on the inspection mission’s objective, the results 


will be tested, according to the following procedures: 
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1. Objective: Reduce the uncertainty of Xj. 


In this case, the decision maker wants to get a more accurate value of the 
fraction of ‘bad” ships in each area. To accomplish that purpose, he/she decides on an 
initial allocation using (2.5) and then continues that allocation. 

After the t™ inspection period is completed, a coefficient of variation, 
CV,(t) , is used as a Statistic to measure the variability or relative dispersion of the 
observed proportion of illegal ships found in each area. CV;(t) is defined as: 

CVj(t) = ott) (2.7) 
L(t) 
Since the proportion of illegal ships after t inspection periods has a beta 


distribution with parameters a(t) and A(t), substitution of Equations (2.1), (2.2) and 


(2.6.a,b) into (2.7) results in the coefficient of variation for the Beta distribution: 


A(t-I+d(t) - k(t) 
[a(t—1) + Bit-1) + d(t) + If [a(t-1) +k] 


Once, the CV;(t) value is under a given limit, g, the decision maker 





CV,(t) = (2.8) 


accepts the updated mean of the proportion of illegal ships for this time period .t, u(t), as 
the actual fraction of illegal ships. How to choose the value of @ is discussed in Chapter 


LI. 


a Objective: Maximize the number of illegal ships found in each area j. 


After observing the results of the t” inspection period the decision maker 
establishes a confidence interval for the odds ratio value € to decide how to reallocate 
assets. Denote the lower and upper confidence limits by LCL and UCL respectively. 


P(X > LCL) 


For any area j if 
P(X < LCL) 


< €, the decision maker will move resources from 


area j to other areas, so as to decrease A; (t+1). 


ie 


P(X > UCL) 


> €, the decision maker will move resources from 
P(X < UCL) 


For any area j if 


other areas to area j, so as to increase A, (tt1). The appropriate value of € to be used is 


explained in detail in Chapter ITI as is the choice of LCL and UCL. 


D. THE ASSESSMENT OF THE PRIOR INFORMATION 


The use of Bayes’ theorem gives an excellent opportunity to exploit all the 
information related to this problem, regardless of its structure or numerical representation. 
However, it is necessary to be careful in the selection of the “prior” intelligence as input 
for the initial assets allocation. First, unless there is strong evidence of the accuracy of the 
intelligence source, the planner should assume a large standard deviation oj(0) for the 
fraction of ships containing illegal cargo in region j. A large standard deviation will allow 
the inspection results in the first period to change the parameters of the beta distribution 
considerably, allowing the estimated value of “4 j;(1) to be much different than the initial 
estimate 4 ;(0). In most, but not all cases, the posterior standard deviation will be smaller 
than the prior standard deviation if many ‘bad” shipments are found. This is intuitively 
reasonable, since the number of ships to be inspected, dj(1), will provide new information, 
and increased information should reduce the uncertainty about the fraction of ships that 
carry illegal shipments, X;. Thus, it might be expected that additional information would 
reduce the standard deviation. 

To illustrate this point, an estimated forecast or intelligence is used to observe the 


sensitivity of the mean 4z;(1) to variations in the prior information. In this example, a 
mean, “£4 ;(0) = 0.16 and a small standard deviation, o ,(0)= 0.06 for X; are assumed, 


which give a(0) =5.8 and £(0) = 30.5 using Equations (2.3) and (2.4). 
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Assume that an inspection of 12 ships is made (1.e., dj(1)=12). The beta posterior 
parameters are calculated assuming various values of the number of inspected ships 


carrying illegal shipments during a time period t, k;(t), and are shown in the Table 1. 





Table 1. Posterior Beta Distribution Parameters Values. Small Variance 

Table 1 shows how the inspection results barely affect the present information, 
given the prior standard deviation. For instance, the difference between the prior mean and 
the posterior mean for eight bad ships found out of 12 ships total is relatively small. 
Therefore the decision maker might conclude that the inspection does not improve the 
information of the proportion of illegal ships. This point confirms the ideas expressed 
above; a prior distribution with a small standard deviation does not lead to a posterior 
distribution that reflects the inspection results. For example, if eight bads were found in 12 


inspections, the raw fraction of bads suggests that the fraction of bads is more likely to be 
closer to of, than 0.28. The eight bads could be a rare event that happened to occur. But 


One would have to be very certain of the initial estimate of 0.16 in order to make this 


interpretation. Figure 3 displays the posterior beta pdf for each number of bad ships found. 


14 


—®— PRIOR 

K valua 0 
—h— K valuc 2 
—>— K valuc 4 


—M— K value 6 


50 pa = 
f(p| dj,k) ; A K value 8 
40 PE 


30 


o vr oOo fF ND rl OoOlUWOlUNNN wo 
ocUCh Cc lOoUlUlCODOWllCOCOClCOllCOFD CCD iC — —_— Q 


Bad Ships proportion (Xj 


5: a — nn: -J 
a ee 
oc So 6 





Figure 3. Posterior Beta pdf distributions. Small Standard deviation 
On the other hand, a prior distribution with a large standard deviation will reflect 
easily any small change in the inspection results. This means that the ‘prior information” 1s 
more sensitive to the inspection results; this indicates the need for inspection and gives to 
the decision maker a better perception of the new information and the value of Xj. 


In this example, a large prior standard deviation o ;(0)= 0.1 is assumed and yz ;(0) 


0.16 as before. The prior distribution with large standard deviation results in an 
unambiguous difference for each number of bad ships found during the inspection of 12 
ships. Table 2 shows the inspection results with their respective posterior beta distribution 
parameters. In this case, there is a substantial difference between the means of the distinct 
inspection results. A larger prior standard deviation results in a easier discrimination of the 


inspection mission results due to the assets allocation. 


| | INSPECTED SHIPS (d, (t))= 
| __|__BETAPOSTERIOR VALUES || 
#BADS(h)/ _aA(1) 


40 | moa | o63 | omnsa7 
60184 | 0.245 | 0.00727 _ 
a 


14.4 0. 409 0. 00950 


eae 





Table 2. Beta Posterior Distribution Parameters values. Large Variance. 
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Figure 4 shows this behavior graphically using the pdf for each Beta Posterior 


distribution values, given the number of bad ships found. 





Figure 4. Posterior Beta Pdf Distributions. Large Standard deviation 

Another useful consideration of Bayes’ theorem is that the decision maker can wait 
and revise the estimated fraction after several inspection periods, using the combined 
inspection results. Given that inspections are independent each other, the same result is 
obtained whether the calculations are made after each inspection period or just once at the 
end of all of the inspection periods. This feature of Bayes’ theorem is very important, since 
it reduces the number of applications of the theorem required to revise the proportion of 
illegal shipments on the basis of various inspections. Of course, the decision maker may 
want to revise after each inspection period, primarily because the resulting proportions 


may help to decide whether or not to execute another inspection. 
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HI. INFERENCE AND MARGINAL ANALYSIS 


A. THE ASSESSMENT OF THE “DECISION” REGION. 


iL Objective A: Reduce the uncertainty of X; 


In Chapter II our first criterion for allocating resources is that the decision maker 
should attempt to estimate the fraction of ships carrying illegal shipments as precisely as 
possible, based upon the coefficient of variation, CV;(t). The coefficient of variation is 
compared to a number @ to evaluate the accuracy of 44(f) , the estimated fraction of ships 
carrying illegal cargo in an area j in a period t. 


0j(t) 
L4(t) 





Our criterion, < @, can be reformulated as a function of the beta distribution 


parameters @ and £. Rearrangement results in : 


o°(t)< p(t). 
Equations (3.1) and (3.2) and substitution result in: 
a(t) G(t) <9 (a(t) 
(a) + AO) (aH+B0+1) — (a+BOY 
Hence, 
A(t) < gp’ @(t)(@(t) + A(t) +1). 
When equality holds, substituting this expression for B,(t) in Equation (2.1) and 
simplifying, results in 


l= (1) - we | 


@(t) = : 1) 
Substituting (3.1) into Equation (2.4) yields: 
1- w(t) - w(t)e’?)(1- uo 
gy = LM = Oe") (I~ we) a 


Q L(t) 
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The a and B values at equality represent the posterior beta parameters for which 


the estimated proportion of illegal ships, j(t), meets the accuracy expressed by 9g, 


according to the coefficient of variation, CVj. Figure 5 1s used to illustrate this concept. 
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Figure 5. Posterior beta parameters as a function of @ 

Figure 5 represents the a and B parameters as a function of the desired accuracy ~ 
of the estimated mean of the proportion of illegal ships in area j, ss(t). In this case a 
44(t) value target of 0.3 was employed. Note that the more accurate the estimated fraction 
of illegal shipments is required (i.e. the smaller the ~ value), the higher o and B values are. 
Therefore, if the decision maker attempts to estimate the fraction of ships carrying illegal 
cargo in an area j with a high precision, then a bigger number of ships should be inspected 
to get a more accurate estimation of sz, .This means a higher number of resources should 
be employed or more inspections should be carried out. 

Following the above discussion, it is not clear that a constant value of @ can be 
used for a large range of s(t) values, or what value of @ to use. Choosing a lower value 
of » will ensure a lower variance on X; relative to s(t), but it will require more 


inspections (and hence time periods) to achieve the desired result. 
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Recall that after t inspection periods X; has a beta distribution with mean sy (t), 
and standard deviation © 44 (t). Figure 6 shows the 95% credible interval for X; for values 
of s(t) between 0.1 and 0.3 when @ is 0.05 and 0.1. These credible intervals are 
calculated using the Excel Version 5 built-in BETAINV function. 
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Figure 6. 95% Credible Interval for 44 as a function of 0. 


Note that the smaller the @ value, the smaller the credible interval. This credible 
interval should be interpreted as follows: there is a probability of .95 that the true fraction 
of ships with illegal cargo in an area J, X;, takes a value between the interval limits. In the 


remainder of this thesis a @ value of 0.1 is used. 


Je Objective B: Maximize the number of illegal incoming ships 
intercepted, k(t). 


P(X) > X) 


In Chapter II, the odds ratio 
is x) 


is suggested as the criterion to decide the 


reallocation of resources among the areas. Let € represent the level by which to measure 


We 


the odds of the fraction of illegal ships, given the inspection results. To used this criterion 
we need to decide on a value of x as well as the value of €. 

Recall that s4(t) plays a central role in determining the allocation of assets using 
Equation (2.5). Therefore we should evaluate how the outcomes of the inspection in any 
period effect the estimation of ,(t). For this reason x will be set equal to sy (t). 

Suppose objective A is complete in period t. Then Xj; has a beta distribution with 
mean 44(t) and standard deviation pyzj(t). The odds that X; exceeds ,zs(t) are 
P(X) = w() 
P(X < (0) 


0.97 for every value of s(t) between 0.02 and 0.3 when o = 0.1. In period (t+1) the 


, and are shown in column 2 in Table 3. It is seen to be between 0.95 and 


number of illegal ships found is denoted by Kj(t+1). When this random variable takes on 


the value kj our odds ratio, based on this new information, 1s PS = wCOIKY = ky) 
P(X; < (OD |Ki(t) = ki) 

For small values of 14 (t) we would not expect large numbers of illegal ships to be 
found. If a large value of k; is found in period (t+1) it would indicate that 4s (t+1) has 
increased from ,4(t), and so more assets should be moved into area j. Similarly, a small 


value of k; in (tt+1) would indicate that 44 (t+1) is smaller than ,4 (t), and so assets can be 
moved out of area j. Columns 3 through 10 in Table 3 show these (posterior) odds for 
values of kj(t+1) between 0 and 7 when dj(t+1) is 12. 

Having set x equal to 44(t), the next problem ts to decide on a value for ¢. As has 
just been pointed out, a high or low conditional odds indicates a change in the value of 
4 (t). For reasons that are presented below we choose an upper value of € to be 1.5 anda 
lower value of 0.67 when s(t). is 0.16. Using these € values in Table 3, if 5 or more 
illegal ships are found in period (t+]1) we conclude that j4(t+1) > s(t) and use a new 


value together with other new values from other areas to reallocate assets using Equation 


(2.5). Similarly, if no illegal ships are found, we conclude that s,(t+l) < s(t) and 


reallocate accordingly. If between | and 4 illegal ships are found we conclude that 4, (t+1) 
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= u(t). The calculation of s;(t+1) when it is considered different from ,4 (t) is shown in 
Section B below. 

To justify the choices of 1.5 and 0.67 for upper and lower values of € when sy (t). 
is 0.16, we further analyze the outcome of the inspection results in period (t+1). At the 
start of period (t+1), it is reasonably to assume that the number of illegal ships found in 
period (t+1) is a binomial random variable with parameters dj(t+1) and 44 (t). For the case 
when 44 (t)= 0.16, P(Ki(t +1) = 5) = 1-B(4; 12, 0.16) = 0.0310 and P(K,(t +1) = 4) = 1- 
B(3; 12, 0.16) = 0.111. Thus using a 5% cut-off value, we would claim 44 (t+1) > s(t) if 
k,(tt1)= 5. This agrees with the value obtained using the € value of 1.5 in Table 3. For 
other values of s4(t) in the range of 0.02 and 0.30, the critical levels of kj(t+1) are 
approximately equal using ¢ values of 1.5 and 0.67 in Table 3, or cut-off values of 0.05 
and 0.95 with the cumulative tail distribution of the binomial. These values are tabulated 
and presented in Appendix B. 

This above example shows the procedure of calculating the required 95% odds 
ratio € bounds. These bounds depend on the s(t) value obtained in the completion of 
objective A and the number of ships to be inspected, dj(t) Given that each area can have its 
own 14(t) value, its own asset inspection capacity and its own rate of incoming ships; the 


odds ratio bounds must be calculated for each area using the procedure described above. 


B. INFERENCE AND DECISION FOR ASSETS ALLOCATION 


To illustrate the preceding ideas, Table 3 is used to present the Odds values 


for different values of 1;(t) from 0.02 to 0.30 and k; from 0 
P( X(t +1) < OIK(t) = ky) 


to 7. Assume objective A) is completed at time t using a equal to 0.1, for some area j 
t(t) = 0.16 , and so oj(t) =0.016. The probability that the true fraction of illegal ships, 
Xj(t), is greater than 1(t) is 0.49 so that the odds are 0.49/0.51= 0.96 (see column 2). 
Suppose that during inspection period (t+1) 6 of the 12 ships inspected are found to be 


carrying illegal goods. Using Equations (2.6.a) and (2.6.b) the new a(t+1) and B\(t+1) are 
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89.974 and 446.865, respectively. Using Equation (2.1) the value of u(t+1) is 0.167, and 
P(X(t +1) > 0.16{Kj(t +1) = 6) =0.674, so the odds are 0.674/0.328 = 2.07, a value 


above the upper bound ¢ of 1.5. The increase in odds from approximately 1 to 2 is a 
strong indication that the true fraction of illegal ships entering area j has increased, and 
that the finding of 6 illegal ships in 12 inspections can not be reasonably explained as 


simply a random occurrence. 
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Table 3. Odds Ratio 





When an unexpected number of ‘bad” ships 1s found in period (t+1), the decision 


maker needs to calculate the s,(t+1) value reflecting the extreme odds ratio value. In 


Chapter II, it was shown that it is necessary to use a quite high variance to detect the 


changes in 44 (t+1) caused by unexpected high or low kj(tt+1) values. Having decided that 
44 (t+1) is different from the jy(t) obtained after completion of objective A, it is 
appropniate to disregard data from periods | through t, and base the estimate of sy (tt+1) 


only on the inspection results found in period (t+1). Using the fraction of ships inspected 


that are found to be illegal we set: 


u(t+1) = a (eaen 


22 


o(t+l) = gu(t+)) (3.3.a) 
where @ takes a value of 0.5, large enough to discriminate the possible changes in 
posterior inspections, based upon the ideas expressed in Section D, Chapter II. Then 
Equation (2.5) is used to estimate the asset allocation to each area according to the new 
Ss. 

Having the updated 44 and oj, the new a; and B; are calculated using Equations 
(2.3) and (2.4). At this point the decision maker restarts the inspection process following 


the 44 estimation procedure under the Bayes procedure and using the conjugate beta 


distribution. 
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IV. MODEL IMPLEMENTATION. 


A. INTRODUCTION 


To test the reasonableness of the model, the decision cycle calculations are 
implemented using an Excel spreadsheet to observe the possible values obtained after an 
inspection is executed, the calculation of the posterior beta parameters and the 
corresponding inference for each X; using the mean value of the posterior beta 
distribution, u;. This representation considers two objectives mentioned above: A) Reduce 
the uncertainty of X;. and B) Maximize the expected number of illegal ships found in each 
inspection. Each objective is tested under four different scenarios. These scenarios 
represent inspection missions of length / of one month. The assets inspection capacity is 
expressed by ships per day. Each scenario is defined by: a) Total number of available 
assets (N), b) True fraction of ships carrying illegal cargo, Xj, c) Accuracy level of the 


estimated mean fraction of ships carrying illegal cargo, g@ value and d) Odds ratio € value. 


Two procedures are considered: 1) The decision maker, using the prior 
information or “intelligence”, decides first to carry out several inspection missions to get a 
more accurate estimate of the fraction of ships with illegal cargo, 4, for each area J. 
Second, he/she allocates assets using the obtained uy, to make inspections with the 
objective of maximizing the total number of ships detected with illegal cargo. 2) The 
decision maker decides not to minimize the uncertainty of the fraction of ships with illegal 
cargo, Wj, aS his initial objective. Instead, the planner chooses to maximize the total 
number of ships detected with illegal cargo, using only the pnor information obtained by 
intelligence. 

This model contains a single random factor in each period for each area: the 
number of illegal ships found in an area j during an inspection period, K,(t). Given this, a 
spreadsheet can be uSed as a simulation tool to demonstrate the analysis and the model 


plausibility. 
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B. EXAMPLE USING PROCEDURE 1. 


Recall that in procedure 1, we start with intelligence estimates as the values of 
14(0) in each area j and first reduce the uncertainty of these values. The model 
representation is displayed in Table 4 for one scenario. The number of inspected ships 
carrying illegal cargo are independent binomial random numbers with probability of 0.2. 
From intelligence reports the fraction of illegal ships in areas 1, 2 and 3 are 0.3, 0.4 and 
0.25, respectively. Following the arguments at the end of Chapter IJ, we choose a 
coefficient of variation equal to 0.5. Using this value the standard deviations for areas 1, 2 


and 3 are 0.15, 0.2 and 0.125, respectively. Other scenarios are tabulated in Appendix A. 


OBJECTIVE: REDUCE UNCERTAINTY OFX; 


Total Assets(n) 25 cp Value Q1 
Real X, §=1,2,3) Q2 


Variable description 
Area information 


Ship arrival rate 

Assets inspection capacity 

Initial mean fraction of illegal ships 

Initial SD. fraction of illegal ships 
Initial Beta parameters 


First Inspection 
Illegal ship rate 
Initial assets allocation 
# of ships to be inspected 
# of illegal ships found 

Updated Beta distribution 
Alpha 
Beta 40.83 
Updated mean fraction of illegal ships Q22 
Updated SD. fraction of illegal ships 007 
Coefficient of variation 258 
Updated illegal rate 197.771 





Table 4. Model Representation. Reduce the uncertainty of X; 
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Second Inspectlon 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. of fraction of illegal ships 
Coefficient of variation 
Third Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D.fraction of illegal ships 
Coefficient of variation 
Fourth Inspection 
New illegal rate 
Assets allocation 


# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 
Fifth Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 


| # of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of Illegal ships 


Coefficient of variation 
Sixth Inspection 


New Illegal rate 
Assets allocation 
# of ships to be inspected 


# of Illegal ships found 
Updated Beta Distributlon 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 


Aj2 
dj2 
K 2 


a,(2) 
B,(2) 
4,(2) 
6,(2) 
CV 


Ajw,(2) 
Aj3 
d,3 

K(3) 


a;(3) 
B,(3) 
(3) 
5,(3) 
CV 


11403) 
A,4 
di4 
K,4 


a,(4) 
B,(4) 
(4) 
o(4) 
Cy 


Ain (4) 
A5 
d\5 
K)5 


a,(5) 
B,(S) 
u,(5) 
3,(5) 
CV 


A,u{5) 
Aé6 
dé 
K6 


a,(6) 
B,(6) 
u,(6) 
6,(6) 
CV 


1 
44 
10 


21.50 
74.83 
0.22 
0.042 
0.189 


198 
1 
44 

7 


28.50 
111.83 
0.2 
0.034 
0.149 


180 
ie 
44 
10 


38.50 
145.83 
0.21 
0.030 
0.143 


183 
1 
ad 


44.50 
183.83 
0.193 
0.0262 
0.134 


175 401 
11 
44 
8 


52.50 
219.83 
0.19 
0.024 
0.124 


17.00 
88.00 
0.16 

0.036 
0.221 


120 
10 
50 
10 


27.00 
128.00 
0.17 
0.030 
0.150 


127.5 
10 
50 

3 


36.00 
163.00 
0.18 
0.027 
0.151 


135 
10 
30 


44.00 
211.00 
O17 
0.0236 
0.137 


129.412 
10 
50 

5 


49.00 
256.00 
0.16 
0.021 
0.131 


Table 4. Model Representation. Reduce the uncertainty of X; (Continuation). 


a] 


16.75 
42.25 
0.28 
0.058 
0.205 


154 
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20.75 
62.25 
0.25 
0.047 
0.178 


137.5 


24 


26.75 
80.25 
0.25 
0.042 
0.167 


137.5 


24 


29.75 
101.25 
0.23 
0.0365 
0.161 


124.905 
4 
24 
4 


33.75 
121.25 
0.22 
0.033 
0.152 


Seventh Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Eigth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of legal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Ninth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Tenth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 
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Ap4(6) 
A? 
d,7 
K,7 


04(7) 
(7) 
Ne 


3,(/) 
CV 


Ag (7) 
As 
d8 


K8 


04(8) 
B;(8) 
}4(8) 
o,(8) 
CV 


p48) 
AQ 
d9 


J 


K9 


04,(9) 
B,(9) 
14(9) 
o,(9) 
CV 


Ap4(9) 
A,10 
d,10 


K,10 


at,(10) 
B,(10) 
p,(10) 
6,(10) 
CV 


61.50 
254.83 
0.19 
0.022 
0.114 


174.974 
11 
44 

7 


68.50 
291.83 
0.19 
0.0206 
0.109 


171.092 
11 
a 

g 


77.50 
326.83 
0.19 
0.020 
0.096 


172.506 
11 
44 

9 


86.50 
361.83 
0.19 
0.019 
0.096 


120.492 


59.00 
296.00 
0.17 
0.020 
0.119 


124.648 


73.00 
332.00 
0.18 
0.0191 
0.106 


135.185 
10 
50 
13 


86.00 
369.00 
0.19 
0.018 
0.092 


141.758 
10 
50 

8 


94.00 
411.00 
0.19 
0.017 
0.093 


119.758 
4 
24 
7 


40.75 
138.25 
0.23 
0.031 
0.137 


125.209 
4 
24 
3 


43.75 
159.25 
0.22 
0.0288 
0.134 


118.534 


24 


47.75 
17825 
0.21 
0.027 
0.121 


115.694 


24 


51.75 
198!25 
0.21 
0.025 
0.124 


Table 4. Model Representation. Reduce the uncertainty of X;_ (Continuation). 


Eleventh Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 


Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Twelfth Inspection 
New illegal rate Agui(1 1) 
Assets allocation A\1 Z 
# of ships to be inspected 1 2 
# of illegal ships found K1 2 

Updated Beta Distribution 
Alpha o, (12) 104.50 110.00 59.75 
Beta B12) 431.83 495.00 239.2 
Updated fraction of illegal ships 14(12) 0.19 O18 020 
Updated SD. fraction of illegal ships o,(12) 0.01709 0.01567 0.02309 
Coefficient of variation CV 0.088 0.086 0.101 





Table 4. Model Representation. Reduce the uncertainty of Xj. (Continuation). 


The following observations are made from the results in Table 4. 


1. CV reduction as a measure of uncertainty of X;. 


Figure 7 illustrates the CV reduction for each area. Notice that for the area with 
bigger initial variance, area 3 , it 1s necessary to carry out more inspection missions to 
achieve a variance that satisfies o, = 0.1u1;. For this example, it took nine months in areas | 
and 2, and 12 months in area 3. This shows the importance of intelligence related to this 
type of operations. The number of time periods required to obtain a good uw, estimation 


depends on how accurate the intelligence is. 


2° 


CV Reduction 


—@CV Area 1 
—8— CV Area 2 
—&— CV Area 3 


——— Limit Value 


CV value 


Time Period (t) 





Figure 7. Coefficient of Variation Reduction 


De All the estimated fraction of ships carrying illegal cargo, p(t) converge 
to the unknown real proportion , 0.2 in all areas. 


Figure 8 displays how the 1,(t)’s values get closer to the real proportion of illegal 


ships in each area This result indicates that the model can be used as an adequate tool to 


get good approximations to the real proportion. 


Mean Value 


Si Oo 7 a8. 9) LO ele 
Time Period (t) 





Figure 8. Estimated fraction of ships carrying illegal cargo, t,,. 


30 


Having reduced the uncertainties on each pj; to within a specified limit, assets are 


now allocated to maximize the number of illegal ships found. 


Table 5 shows the process using this objective. Equation (2.5) is used to 


recalculate the assets allocated to each area j at the start of every period. 


OBJECTIVE: MAXIMIZE NUMBER OF ILLEGAL SHIPS FOUND 


Total Assets (n) 
Real X, Q=1,2,3) 
Variable description 


Area information 

Ship arrival rate 

Assets inspection capacity 

Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 


Initial Beta parameters 
Alpha 


Beta 

First Inspection 

[legal ship rate 

Initial assets allocation 

# of ships to be inspected 

Expected number of illagal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 

Updated illegal rate 

Second Inspection 

New assets allocation 

# of ships to be inspected 

Expected number of ilegal ships 

# of illegal ships found 





95% ODDS RATIO BOUNDS Area 1 Area 2 Area 3 


LOWER 0.36 0.43 0.44 
90 1.99 1.96 1.72 


Symbol Area 1 Area 2 


A, 900 750 
Cap, 4 5 


},(0) 
6,(0) 


0,(0) 
BO) 


Agu(0) 
A{1) 

dj(1) 

E[Kj(1)) 

K(1) 
P(X>p,(0))/PCX<4,(0)) 
},(1) 

o,(1) 


a(1) 
BC) 
Ayyi(1) 


A2 
dz 


K2 


Table 5. Maximizing the number of illegal ships found. 


at 


Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


Third Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 

New illegal rate 


Fourth Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Update Beta Distribution 
Alpha 

Beta 

New illegal rate 


Fifth Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Update Beta Distribution 

Alpha 

Beta 


New illegal rate 


PX>y,(1))PCX<p(1)) 
y,(2) 
6,(2) 


a,(2) 
B;(2) 
Aju4(2) 


A3 
d3 
E[Kj(3)} 

K,(3) 
P(X>H(2))/P(X}<p,(2)) 
H,(3) 

6,(3) 


o4(3) 
B,(3) 
AHy{3) 


A,4) 
d,4) 
E{Ky(4)] 

K4) 
PX>H,(3))/PCX<H,(3)) 
14,(4) 

0,(4) 


a,,(4) 
B,(4) 
Ayui(4) 


A,9) 
d\3) 
ETKj(5)] 

Ky9) 
P(X>14,(4)) PCy<u,(4)) 
H,(5) 

6,(5) 


04(5) 
BO) 
Ajuj(5) 


0.94 
0.19 
0.019 


80.81 
344.51 
171 


12 
48 
9.6 


0.65 
0.19 
0.019 


80.81 
344.51 
171 


48 
9.6 


0.44 


0.19 
0.02 


80.81 
344.51 
171 


1 

44 

8.8 
8 


0.90 
0.19 
0.02 


80.81 
344.51 
171 


112 
0.18 
0.018 


81.82 
372.74 
135 


40 
8.0 


0.93 
0.18 
0.018 


81.82 
372.74 
135 


40 
8.0 


0.93 


0.18 
0.02 


81.82 
372.74 
135 


39 
7.0 


0.91 
0.18 
0.02 


81.82 
372.74 
135 


0.79 
0.20 
0.019 


79.80 
319.20 
110 


30 
6.0 


0.44 
0.20 
0.020 


79.80 
319.20 
110 


30 
6.0 

10 
2.09 


0.33 
0.17 


2.33 
4.67 


183 


8 
48 
19.2 
13 
0.22 


0.27 
0.14 


2.65 
12 
149 


Table 5. Maximizing the number of illegal ships found. (Continuation) 


a2 


The observations obtained from this simulated decision process are the following: 


3. Observation 3. The odds ratio criterion 


Table 6 presents the odds ratio value obtained for the example presented in Table 
5. Except for area 3 during the fourth and fifth inspection, all the odds values are inside 
the ¢ values bounds,. From the results, it is clear that when the odds ratio value is out of 
bounds the model calculates the new assets allocation as a response of the change in the 
estimated fraction of ships carrying illegal cargo, p;. Even though a small change might be 
considered “normal” due to the process randomness, the model estimation is quite robust 
and calculates the appropniate asset allocation as result of the new 1;’s values. Any odds 
ratio value out of limits might be used as an indicator or warning for assets allocation in 
future inspection periods. In fact, the decision maker may select a desired X; value and use 


the odds ratio as a “sensor” value to detect a change of unlawful events of this nature. 


POSTERIOR ODDS RATIO 
INSPECTION _—_ Areal ‘Area2 Area 3 
1ST 0.645 0.926 0.793 
2ND 0.941 1.116 0.793 


3RD 0.65 0.93 0.44 
4TH 0.438 0.926 2.087 
STH 0.90 0.91 0.22 





Table 6.Posterior Odds Ratio. 


4, Observation 4: Maximum number of illegal ships found. 


Table 7 displays the number of ships found during the 5 inspection periods for the 
example in Table 5. Note that the number of illegal ships found 1s close to the predicted 
expected number of illegal ships calculated according to the updated 1, values. These 
results are quite satisfactory under the assumptions presented in this thesis. They 


demonstrate the feasibility of getting acceptable outcomes given the number of assets 


available. 
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EXPECTED NUMBER OF ILLEGAL SHIPS 
INSPECTION Area 1 Area 2 


IST 7 
2ND 
3RD 
4TH 
STH 


# ILLEGAL SHIPS FOUND 
INSPECTION Area 1 Area 2 Area 3 Total 


IST 7 7 5 


2ND 9 
3RD 7 
4TH 5 
STH 8 





Table 7.Expected Number of Illegal ships vs Number of Illegal ships found. 


C. EXAMPLE USING PROCEDURE 2. 


Sometimes the decision maker is constrained by time and other operational 
requirements. So, he/she decides to execute inspection missions to enforce laws trying to 
capture the maximum number of ships with illegal cargo starting with the intelligence of 
the j1;’s in each area. Recall from Table 4 that our initial estimates are 0.3, 0.4 and 0.25 for 
areas 1, 2 and 3 respectively. Table 8 presents the results of a simulated inspection process 
where the estimates of aj(t), B,(t), p(t) and oj(t) are obtained using the Bayesian 


procedure. The assets are reallocated using Equation (2.5). 
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OBJECTIVE: MAXIMIZE NUMBER OF ILLEGAL SHIPS FOUND 





Total Assets (n) 
Real X, 


Variable description 


Area information 


Ship arrival rate 
Assets inspection capacity 


Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 
Initial Beta parameters 

Alpha 

Beta 

First Inspection. 

Illegal ship rate 

Initial assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Updated Beta parameters 

Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated illegal rate 


Second Inspection 

New assets allocation 

# of ships to be inspected 
Expected number of illegal ships 
# of illegal ships found 


Updated Beta parameters 
Alpha 


Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


New illegal rate 


Third Inspection 


New assets allocation 
# of ships to be inspected 
Expected value of illegal ships 


# of illegal ships found 


Updated Beta parameters 
Alpha 


Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


New illegal rate 


Symbol 


hj 
Cap, 
H,(0) 
o,(0) 


at,(0) 
B,(9) 


AH {0) 
At) 
d,(1) 

E(Kj(1)] 
K(1) 


aj(1) 
B,(1) 
H(1) 
o,(1) 
Ap (1) 


A,2) 
d,(2) 
E{Kj(2)] 
K 2) 


at,(2) 
B,(2) 
p1,(2) 
62) 
hp {2) 


A,3) 
d,3) 
E(K\(3)] 
K(3) 


a,(3) 
B(3) 
(3) 
o(3) 
AH (3) 


Areal Area2 Area3 


900 
4 


0.30 
0.1500 


2.50 
5.83 


270 
in 
44 
8.8 

9 


11.50 
40.83 
0.22 
0.06 


197.77 


20.50 

75.83 
0.21 
0.04 

181.52 


30.50 
109.83 
0.22 
0.03 
195.61 


750 
5 


0.40 
0.2000 


2.00 
3.00 


300 
10 
50 
10.0 
12 


14.00 
41.00 
0.25 
0.08 


190.91 


45 
9.0 
12 


26.00 
74.00 
0.26 
0.04 


195.00 


9 
45 
9.0 
14 


40.00 
105.00 
0.28 
0.04 
208.80 


550 
6 


0.25 
0.1250 


2.15 
8.25 


138 
4 
24 
4.8 
5 


7.75 
27.25 
0.22 
0.07 


121.79 


30 
6.0 


13.75 
51.25 
0.21 
0.05 

116.35 


30 
6.0 


18.75 
76.25 
0.20 
0.04 
108.55 


Table 8. Maximizing the number of illegal ships found, using only “intelligence”. 


Fourth Inspection 

New assets allocation 

# of ships to be inspected 
Expected value of illegal ships 
# of illegal ships found 


Update Beta Distribution 

Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


New tlegal rate 


Fifth Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Update Beta Distribution 

Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


New illegal rate 


Table 8.(Continuation). 


A,4) 
d4) 
E(Kj(4)} 
K\4) 


a,(4) 
B,(4) 
H,(4) 
o,(4) 
Au {4) 


AyS) 
dj9) 
E[Kj(5)) 
KS) 


a,(3) 
B,(S) 
H,(5) 
6,(5) 
Aiu,(5) 


10 


50.00 
145.00 
0.28 
0.03 
192.31 


58.00 26.75 
184.00 118.25 
0.23 0.19 
0.03 0.03 
175.00 102.88 





From this examples the remarks obtained are the following: 


Ie With high uncertainty of p;., the odds criterion can not be used for 
asset allocation. 


Table 9 displays the asset allocation and the mean of the fraction of illegal ships 
carrying illegal cargo from the Bayesian procedure. The initial high vanability or 
uncertainty of Xj produces a very wide confidence interval. Therefore, after inspection 
results, the obtained odds value does not offer clear information for asset allocation. If it is 
used would cause large changes in the asset allocation. 


ASSETS 


Area 2 


MEAN 


INSPECTION Area 1 Area 3 INSPECTION Area 1 Area 2 Area 3 


IST 


2ND 


3RD 
4TH 


ar 


1 
13 
1 
10 
i's 


10 
g 

10 
1 
10 


IST 
2ND 
3RD 
4TH 


Sie 


Table 9. Assets Allocation. 
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0.239 
0.168 
0.165 
0.183 
0.192 


0.236 
0.220 
0.260 
0.234 
0.220 


0.136 
0.146 


0.166 
0.146 


0.150 





Even though the odds value criterion is not appropiate for this situation, the 
decision maker can decide the asset allocation by the weighted distribution expressed by 
Equation (2.5). In this case, the algorithm offers a good alternative when the planner 
desires to get a satisfactory inspection result and he/she can not obtain a good estimate of 
the fraction of ships that are transporting illegal cargo. Table 10 shows a comparison of 
the expected number of illegal ships and the number of illegal ships found for the example 


in Table 8. The results are quite acceptable, in spite of the inaccurate initial estimation of 


Lj. 


EXPECTED NUMBER OF ILLEGAL SHIPS # ILLEGAL SHIPS FOUND 
INSPECTION Area 1 Area 2 Area3 Total Areal Area2 Area3 Total 
IST 9 5 12 
2ND 


3RD 
4TH 
STH 





Table 10. Expected number of illegal ships and number of illegal ships found. 


Pe. Inspection results guide to the real value of X; 


Figure 9 displays the estimated 1 ; values for each of the three regions. Even 
though the planner decides to carry out inspections with a rough estimation of uj, he can 


obtain estimated 1; values with tendency to the real value of Xj. 


— MEAN Area 1 —— MEAN Area 2 —@— MEAN Area 3 MEAN Real Xj 


0.40 & 
0.35 
0.30 
0.25 
0.20 
OD 
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0.00 
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= 
a 
- 
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| 
w 
— 
&, 
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Figure 9.Estimated value of the fraction of ships with illegal cargo. 


oy 


The asset allocation variability is low through the different inspection periods. 

A comparison between the courses of action that a decision maker can execute 
using this model, drives the important observation that despite the decision maker uses 
Equation (2.5) under unreliable 1;’s values, the assets allocation keeps a regular pattern in 
each area during the time periods. Table 11 illustrates this point. After the first inspection, 
only the assets corresponding to area 2 present a significant change in the asset allocation. 
However, there exist an “equilibrium” after the inspections results update the estimate of 


the ,1;’s values. 


ASSETS 


INSPECTION Area 1 Area 2 Area 3 
1ST 11.0 10.0 4.0 


2ND 11.0 9.0 9.0 


3RD 11.0 9.0 9.0 
4TH 11.0 10.0 4.0 


STH 11.0 9.0 4.0 





Table 11. Assets allocation. 


3. Observation 4: Case 1 and Case 2 trade-offs. 


In Case 2, unlike Case 1, the decision maker consider reallocation of assets 
between different assets of the end of every period t. In case 1, the asset allocation is fixed 
when the decision maker desires to get precise u1;'s values for each area. Then he/she 
reallocates assets to obtain the maximum number of ships carrying illegal cargo. Recalling 
that the underlying goal of this type of operations 1s to capture those ships involved in this 
kind of unlawful activity, it is useful to see what course of action offers the best tradeoff, 
considering time and the uncertainty involved in the value of the fraction of ships 
transporting illegal shipments, X;. To analyze this point, Table 12 is used to illustrate the 
number of illegal ships found using both procedures. Note that the average number of 
illegal ships found per inspection period is slightly higher using procedure 2. From a 


practical point of view, the results show that using procedure 2 produces at least as good 
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results as procedure 1. 111 illegal ships were found in the first five periods using 


procedure 1, whereas this number was 123 using procedure 2. 


Number of | Number of | Number of Illegal 





Table 12. Number of Illegal ships found through Procedure 1 and Procedure 2. 

In conclusion, it is reasonable to conclude that procedure 2 is better from a 
practical consideration. However, procedure 1 represents a useful method to gather 
important information or intelligence related with contraband or illegal shipment. It gives a 
clear insight about the presence of these activities. As it was expressed before, an accurate 
estimate of the fraction of ships transporting illegal cargo gives the decision maker a 
favorable starting point to ensure optimal results in the achievement of Control Maritime 


missions. 


a 





V. CONCLUSIONS AND FURTHER RESEARCH. 
A. CONCLUSIONS 


The model developed in this thesis provides an analytical procedure to determine 
the best allocation of assets for the execution of a Maritime Control mission in closed 
areas. Asset allocation is based on the fraction of ships that are carrying illegal cargo in 
each area j during a time period t. This model uses forecast information or “intelligence” 
to deploy such assets according to a Bayesian estimation of the illegal fraction. Few 
quantitative data related to this issue are available, which makes difficult the use of 
classical statistical methods. 

The model estimates the proportion or fraction of ships carrying illegal cargo 
which heads to harbor through natural straits or channels. Using the available 
‘{ntelligence” as input to generate an initial estimate of this fraction, the model defines a 
decision making process as a guide to the planner for asset allocation. The initial allocation 
decision is based upon a weighted distribution using the initial estimated mean of the 
fraction of illegal ships with illegal cargo, pj. Those areas to be inspected with higher 
fractions will receive more assets. The model demonstrates that if the decision maker 
decides to get a precise estimation of uy; by carrying out several inspection with a fixed 
number of assets in each area, then he/she can maximize the number of illegal ships found 
with little change in the asset allocation. On the contrary, when the planner decides to 
make inspections only with the prior information from intelligence, he/she may reallocate 
assets more frequently, and more inspections are needed to get a desired precision for the 
estimated 4; values. Even though this thesis keeps separate the logistics considerations 
related to the assets allocation, the fact remains that reducing the uncertainty of the p, 
value should be helpful in minimizing the complexity of the allocation problem. 

This thesis involves intelligence as a quantitative input into the decision making 
process. The Bayesian approach makes possible the use of “Subjective” information into 
the analytical model. The model is used to demonstrate that ‘intelligence’ can be 


translated as an important estimator in this stochastic analysis. 


4] 


In summary, this model represents a satisfactory starting point to obtain insights 
into the problem here presented. It presents an adequate alternative to improve the 


planning process of the type of Maritime Control missions presented. 


B. FURTHER RESEARCH 


Many of the model assumptions represent interesting areas for further 


investigations and potential model enrichment. For example: 


I Time to inspect a ship 


The formulation of a stochastic model could provide an excellent tool to estimate 
this parameter. The associated randomness with it makes difficult the determination of the 
required number of assets for an inspection mission. Sea state, ship type to be inspected, 


assets capability, training, etc are some of the factors involved in this issue. 


2 Probability of transporting illegal cargo 


For this study, this value was assumed constant for every ship. However, it is 
known that some ships are more likely than others for transporting illegal cargo. Also, 
preceding data or frequency of ship type involved in this type of outlaw activities, may be 
used to develop a forecast model for determining this probability so that ships can be 


ranked based on this value. 


3. Incorporation of Logistics restrictions 


The scheduling or assets allocation implies constraints or considerations based on 
maintenance, personnel rotation, budgets and other conditions which influence ship 
distribution. Therefore, optimization models are potential tools to evaluate these factors 


and to give additional inputs to the planning of Maritime Control missions. 
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4. Iilegal shipment detection 


The design of a stochastic model or an optimization model might be possible to 
estimate the probability of detecting an illegal shipments on board. The measurement of 
this value represent an important enhancement to the decision making process related to 


Maritime Control planning. 
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APPENDIX A. MODEL IMPLEMENTATION SCENARIOS 


In this appendix the additional model implementation scenarios are 
collected to show how the model evaluates fictitious Maritime Control missions. 

The scenarios are defined by the following factors: 

a) Number of Assets availables, N. 

b) Real fraction of illegal ships transporting illegal cargo in each 
area, X;. However, in this example, it was assumed that this value is equal for all the areas. 

c) Accuracy level of the estimated mean of ships carrying illegal 
cargo, lj. 

d) Odds ratio confidence interval, ¢ value. 

Also, after the completion of each mission, the results are tabulated 


to illustrate the inference of results and other considerations. 
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OBJECTIVE: REDUCE UNCERTAINTY OF X ; 


Total Assets(n) 10 
Real (ab 2) 0.2 


Variable description 
Area information 


Ship arrival rate 
Assets inspection capacity 
Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 
Initial Beta parameters 
Alpha 
Beta 
First Inspection 
Illegal ship rate 
Initial assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 
Updated illegal rate 
Second Inspection 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. of fraction of illegal ships 
Coefficient of variation 
Third Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 


~ Value 


0.1 


Symbol 


Aj 
Cap, 
H;(0) 
G;(0) 


a;(0) 
B;(0) 


A,u(0) 
A 


1 
Kj(1) 


aj(1) 
B,]) 
H(1) 
o,(1) 
CV, 
Aju) 


A2 
d2 
K;2 


a(2) 
B;(2) 
H4(2) 
6;(2) 
CV 


Aw (2) 
A;3 
d3 

Ki(3) 


Table 13.Procedure 1. Scenario 1 
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Areal Area2 


900 780 
4 5 
0.3 0.4 
0.010 0.015 


629.70 426.27 
1469.30 639.40 


270 300 
4 4 
16 20 
0 3 


629.70 429.27 
1485.30 656.40 


0.30 0.40 
0.010 0.015 
0.033 0.038 

267.957 297 
4 ~ 
16 20 
6 3 


635.70 432.27 
1495.30 673.40 


03 0.39 

0.010 0.015 

0.033 0.038 

270 ~=—«-293 
4 4 
16 20 
3 g 


Area 3 


550 


0.25 
0.020 


116.94 
350.81 


119.94 
359.81 
0.25 
0.020 
0.079 


138 


122.94 
368.81 
0.25 
0.020 
0.078 


138 


12 


Updated Beta distribution 


Table 13.Procedure 1. Scenario | 
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Alpha aj(3) 638.70 437.27 124.94 

Beta Bi(3) 1508.30 688.40 378.81 

Updated mean fraction of illegal ships (3) 0.3 0.39 0.25 

Updated S.D.fraction of illegal ships G;(3) 0.010 0.015 0.019 

Coefficient of variation CV 0.039 0.047 0.088 
Fourth Inspection 

New illegal rate Au(3) 270 292.5 137.5 

Assets allocation Aj4 4 A 2 

# of ships to be inspected dj4 16 20 12 

# of illegal ships found Kj4 4 4 9 
Updated Beta Distribution 

Alpha at;(4) 642.70 441.27 126.94 

Beta Bi(4) 1520.30 704.40 388.81 

Updated mean fraction of illegal ships 11;(4) 0.3 0.39 0.25 

Updated S.D. fraction of illegal ships G;(4) 0.010 0.014 0.019 

Coefficient of variation CV 0.033 0.037 0.077 
Fifth Inspection 

New illegal rate Aiu,(4) 270 §=.292.5 =: 137.5 

Assets allocation Ajo 4 4 2 

# of ships to be inspected djs 16 20 12 

# of illegal ships found Kj 6 5 1 
Updated Beta Distribution 

Alpha atj(5) 648.70 446.27 127.94 

Beta Bi(5) 1530.30 719.40 399.81 

Updated mean fraction of illegal ships 1(5) 0.30 0.38 0.24 

Updated S.D. fraction of illegal ships 6;(5) 0.0098 0.0142 0.0186 

Coefficient of variation CV 0.033 0.037 0.077 
Sixth Inspection 

New illegal rate Ayu(5) 267.935 287.132 133.331 

Assets allocation A6 & 4 2 

# of ships to be inspected di6 16 20 12 

# of illegal ships found K6 3 3 2 
Updated Beta Distribution 

Alpha at;(6) 651.70 449.27 129.94 

Beta Bi(6) 1543.30 736.40 409.81 

Updated mean fraction of illegal ships 11,(6) 0.30 0.38 0.24 

Updated S.D. fraction of illegal ships 6;(6) 0.010 0.014 0.018 

Coefficient of variation CV 0.033 0.037 0.076 


Seventh Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Figth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Ninth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illcgal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updatcd mean fraction of illegal ships 
Updated S.D. fraction of illcgal ships 


Coefficient of variation 
Tenth Inspection 


New illcgal rate 
Assets allocation 
# of ships to be inspectcd 


# of illegal ships found 
Updated Beta Distribution 


Aip{6) 
Aj 
47 
K7 


a;(7) 
BY) 
H(7) 
G;(7) 
Cy 


Au{7) 
Aj8 
d8 
K8 


a(8) 
Bi(8) 
H(8) 
o;(8) 
Cy, 


Aw(8) 
AQ 
d9 
K\9 


a;(9) 
Bi(9) 
H(9) 
0,(9) 
CY 


Aiu,(9) 
A\10 
d;10 
K;10 
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267.212 
4 
16 
4 


655.70 
1555.30 
0.30 
0.010 
0.033 


266.906 
4 
16 
4 


659.70 
1567.30 
0.30 
0.0097 
0.033 


266.605 
4 
16 
1 


660.70 
1582.30 
0.29 
0.010 
0.033 


265.105 
4 
16 
1 


284.186 
4 
20 
s 


452.27 
753.40 
0.38 
0.014 
0.037 


281.338 
4 
20 
4 


456.27 
769.40 
0.37 
0.0138 
0.037 


279.195 
4 
20 
4 


460.27 
785.40 
0.37 
0.014 
0.037 


277.121 
10 
90 

7 


132.405 
2 
12 
0 


129.94 
421.81 
0.24 
0.018 
0.077 


129.525 
Z 
12 
4 


133.94 
429.81 
0.24 
0.0179 
0.075 


130.671 
2 
2 
3 


136.94 
438.81 
0.24 
0.018 
0.074 


130.813 
2 
12 
0 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Eleventh Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Twelfth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


a,(10) 
B;(10) 
14,(10) 
o;(10) 
CV 


Aju,(10) 


A\11 
41 
K11 


o(1 1) 
Bd.) 
H(11) 
ca 
CV 


Ajp,(11) 


A12 
di12 
K12 


(12) 
Bj(12) 
H(12) 
o,(12) 
CV 
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661.70 
1597.30 
0.29 
0.010 
0.033 


263.625 
4 
16 
1 


662.70 
1612.30 
0.29 
0.010 
0.033 


262.167 


666.70 
1624.30 
0.29 
0.00949 
0.033 


467.27 
828.40 
0.36 
0.013 
0.037 


270.479 
10 
50 
10 


477.27 
868.40 
0.35 
0.013 
0.037 


266.002 
10 
50 
18 


495.27 
900.40 
0.35 
0.0128 
0.036 


136.94 
450.81 
0.23 
0.017 
0.075 


128.142 
2 
12 
4 


140.94 
458.81 
0.23 
0.017 
0.101 


129.247 


143.94 
467.81 
0.24 
0.01714 
0.101 


OBJECTIVE: MAXIMIZE NUMBER OF ILLEGAL SHIPS FOUND 


Total Assets (n) 
Real X; (j=1,2,3) 
Variable description 
Area information 


Ship arrival rate 

Assets inspection capacity 

Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 


Initial Beta parameters 

Alpha 

Beta 

First Inspection 

Illegal ship rate 

Initial assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 

Updated illegal rate 

Second Inspection 

New assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


10 
0.2 


95% ODDS RATIO BOUNDS 
LOWER 
85 


Symbol! 


Aj 
Cap, 
H4;(0) 
o;(0) 


a;(0) 
B;(0) 


Aip{0) 
A(t) 

dj(1) 

E{Kj(1)] 

K\(1) 
P(X;>p,(0))/P(Xj<1;(0)) 
H(1) 

G,(1) 


aj(1) 
Bi) 
Au(1) 


A2 
d2 
E{Kj(2)] 

K,2 
P(X}>p,(1))/P(Xj<H,(1)) 
H,(2) 

o,(2) 


a;(2) 
B;(2) 
Aipi(2) 


Area l 


0.69 
4.98 


Area I 


900 


4 
0.28 
0.028 


Meg 2 
184.42 


252 


20 
4.0 


1.31 
0.28 
0.03 


71.72 
184.42 
252 


20 
4.0 


2.01 
0.28 
0.028 


71.72 
184.42 
202 
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BZ 


Area 2 


0.40 
2.60 


Area 2 


750 


5 
0.34 
0.034 


65.66 
127.46 


250 


20 
4.0 


1.02 
0.34 
0.034 


65.66 
127.46 
250 


20 
4.0 


129 
0.34 
0.034 


65.66 
127.46 
Z00 


Area 3 


0.83 
3.52 


Area 3 


550 


6 
0.23 
0.023 


76.77 
257.01 


127 


12 
2.4 


1.52 
0.23 
0.023 


76.77 
257.01 
127 


be 
2.4 


1.52 
0.23 
0.019 


76.77 
257.01 
127 


Third Inspection 


New assets allocation 

# of ships to be inspeeted 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraetion of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


Fourth Inspection 


New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraetion of illegal ships 
Updated S.D. fraetion of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


Fifth Inspection 


New assets alloeation 

# of ships to be inspeeted 
Expected value of illegal ships 
# of illegal ships found 
Posterior Odds Ratio 


Updated mean fraetion of illegal ships 
Updated S.D. fraetion of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


Aj3 5 
dj3 20 
E[Kj(3)] 4.0 
K(3) 7 
P(X>14(2)/P(Xi<H,(2)) 1.31 
1;(3) 0.28 
o;(3) 0.028 
a(3) aMiz 
B;(3) 184.42 
Aju (3) 252 
A;4) 5 
dj4) 20 
E{Kj(4)] 4.0 
K\4) 5 
P(X>H,(3))/P(X)<n,(3)) 0.85 
1;(4) 0.28 
5|(4) 0.03 
a;(4) 7AeiZ 
Bi(4) 184.42 
A,p(4) 252 
Aj5) 3 
di.5) 12 
E[Kj(5)] 2.4 
Ko) 3 
P(X,>p,(4))/P(Xj<p,(4)) 0.90 
H,(5) 0.28 
G;(5) 0.03 
a;(5) 71.72 
Bi(5) 184.42 
Au{5) 252 
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20 
4.0 


1.02 
0.34 
0.034 


65.66 
127.46 


255 


20 
4.0 


1.02 


0.34 
0.03 


65.66 
127.46 
255 


3.0 


0.75 
0.34 
0.03 


65.66 
127.46 
255 


12 
2.4 


0.83 
0.23 
0.023 


76.77 
257.01 


127 


2 

2.4 

10 
4.39 


0.83 
0.42 


76.77 
257.01 
458 


24 

9.6 

13 
0.00 
0.54 
0.27 


1.29 
1.09 
298 


ait = 8=«gomp 


go oa vy OO oO BS WL PS 


andes 
SO = & 


ao © ~4N OO Of WO HS — O&O -— 


and 0 meh 
nS = & 


CV Area 1 
0.033 
0.033 
0.039 
0.033 
0.033 
0.033 
0.033 
0.033 


0.033 
0.033 


0.033 
0.033 


Mean Area 1 
0.3 
0.30 
0.3 
0.3 


0.3 
0.30 


0.30 
0.30 
0.30 
0.29 
0.29 
0.29 
0.29 


CV REDUCTION 


CV Area 2 


0.038 
0.038 
0.047 
0.037 
0.037 
0.037 
0.037 
0.037 


0.037 
0.037 


0.037 
0.036 


MEAN 


Mean Area 2 


0.4 
0.40 
0.39 
0.39 


0.39 
0.38 


0.38 
0.38 
0.37 
0.37 
0.36 
0.35 
0.29 


CV Area 3 
0.079 
0.078 
0.088 
0.077 
0.077 
0.076 
0.077 
0.075 


0.074 
0.075 


0.101 
0.101 


Mean Area 3 
0.25 
0.25 
0.25 
0.25 


0.25 
0.24 


0.24 
0.24 
0.24 
0.24 
0.23 
0.23 
0.24 
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INSPECTIONS MISSIONS RESULTS 


INSPECTION 
IST 
2ND 
3RD 
4TH 
STH 


INSPECTION 
IST 
2ND 
3RD 
4TH 
Ut 


INSPECTION 
IST 
2ND 
3RD 
4TH 
STH 


INSPECTION 
IST 
2ND 
3RD 
4TH 
5TH 


INSPECTION 
IST 
2ND 
3RD 
4TH 
STH 


Area 1 
5.0 
5.0 
5.0 
5.0 
3.0 


Area 1 
1.308 
2.014 

1.31 
0.852 
0.90 


Area I 
0.280 
0.280 
0.280 
0.280 
0.280 


EXPECTED 


Area 1 


7 


So ff Lf CO 


Area 1 


7 


wo oO ~s © 


SCENARIO: 
ASSETS 
Area 2 Area 3 
4.0 2.0 
4.0 2.0 
4.0 2.0 
4.0 2.0 
3.0 40 
POSTERIOR ODDS RATIO 
Area 2 Area 3 
1.022 523 
1.288 1.523 
1.02 0.83 
1.022 4.389 
0:75 0.00 
MEAN 
Area 2 Area 3 
0.340 0.230 
0.340 0.230 
0.340 0.230 
0.340 0.833 
0.340 0.542 
NUMBER OF ILLEGAL SHIPS 
Area 2 Area 3 
7 5 
8 5 
4 2 
4 2 
3 10 


# ILLEGAL SHIPS FOUND 


Area 2 


7 


-_ ny Ss © 


Area 3 
9 
5 
2 
10 
13 
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OBJECTIVE; REDUCE UNCERTAINTY OF X ; 


Total Assets(n) 15 (p Value 0.1 
Real X; G=1,2,3) 0.2 

Variable description Symbol Areal Area2 Area3 
Area information 

Ship arrival rate a 900 750 550 

Assets inspection capacity Cap, 4 5 6 

Initial mean fraction of illegal ships 14;(0) 0.3 0.4 0.25 

Initial S.D. fraction of illegal ships G;(0) 0.010 0.015 0.020 
Initial Beta parameters 

Alpha at;(0) 629.70 426.27 116.94 

Beta B;(0) 1469.30 639.40 350.81 
First Inspection 

Illegal ship rate A,u.(0) 270 300 137.5 

Initial assets allocation Al 7 6 2 

# of ships to be inspected dj1 28 30 12 

# of illegal ships found K,(1) 4 12 2 
Updated Beta distribution 

Alpha o(1) 633.70 438.27 118.94 

Beta Bi(1) 1493.30 657.40 360.81 

Updated mean fraction of illegal ships W(1) 0.30 0.40 0.25 

Updated S.D. fraction of illegal ships o,(1) 0.010 0.015 0.020 

Coefficient of variation CV 0.033 0.037 0.079 

Updated illegal rate A(t) 268.138 300 136 
Second Inspection 

Assets allocation Aj2 7 6 2 

# of ships to be inspected dj2 28 30 12 

# of illegal ships found K,2 5 6 3 
Updated Beta Distribution 

Alpha at;(2) 638.70 444.27 121.94 

Beta B,(2) 1516.30 681.40 369.81 

Updated mean fraction of illegal ships (2) 0.3 0.39 0.25 

Updated S.D. of fraction of illegal ships 6;(2) 0.010 0.015 0.019 

Coefficient of variation CV 0.033 0.037 0.078 
Third Inspection 

New illegal rate Ayu,{2) 270 = -293 138 

Assets allocation Aj3 7 6 2 

# of ships to be inspected dj3 28 30 12 

# of illegal ships found K((3) 4 6 4 
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Updated Beta distribution 
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a 


Alpha a;(3) 642.70 450.27 125.94 

Beta Bi(3) 1540.30 705.40 377.81 

Updated mean fraction of illegal ships 11;(3) 0.29 0.39 0.25 

Updated S.D.fraction of illegal ships 6;(3) 0.010 0.014 0.019 

Coefficient of variation CV 0.039 0.046 0.088 
Fourth Inspection 

New illegal rate Aju{3) 261 292.5 = 137.5 

Assets allocation Aj4 7 6 2 

# of ships to be inspected dj4 28 30 12 

# of illegal ships found K;4 4 4 3 
Updated Beta Distribution 

Alpha a;(4) 646.70 454.27 128.94 

Beta Bi(4) 1564.30 731.40 386.81 

Updated mean fraction of illegal ships 11;(4) 0.29 0.38 0.25 

Updated S.D. fraction of illegal ships G;(4) 0.010 0.014 0.019 

Coefficient of variation CV 0.033 0.037 0.076 
Fifth Inspection 

New illegal rate Aju, (4) 261 285 137.5 

Assets allocation Ajo 7 6 2 

# of ships to be inspected djs 28 30 12 

# of illegal ships found Kj9 4 6 3 
Updated Beta Distribution 

Alpha o,(5) 650.70 460.27 131.94 

Beta Bi(5) 1588.30 755.40 395.81 

Updated mean fraction of illegal ships L14;(5) 0.29 0.38 0.25 

Updated S.D. fraction of illegal ships G;(5) 0.0096 0.0139 0.0188 

Coefficient of variation CV 0.033 0.037 0.075 
Sixth Inspection 

New illegal rate Aw(5) 261.559 283.959 137.5 

Assets allocation Aj6 7 6 2 

# of ships to be inspected di6 28 30 12 

# of illegal ships found K6 4 3 2 
Updated Beta Distribution 

Alpha o;(6) 654.70 463.27 133.94 

Beta Bi(6) 1612.30 782.40 405.81 

Updated mean fraction of illegal ships 14,(6) 0.29 O37 0.25 

Updated S.D. fraction of illegal ships 6;(6) 0.010 0.014 0.019 

Coefficient of variation Cy 0.033 0.037 0.075 


Seventh Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Eigth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Ninth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Tenth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 


d:7 
K;7 


o;(7) 
B;(7) 
H(7) 
6;(7) 
CV 


Au(7) 
Aj8 
d8 
K,8 


a;(8) 
B;(8) 
H4(8) 
o,(8) 
Cv 


AjH(8) 
AQ 
dg 
K;9 


aj(9) 
B;(9) 
H(9) 
6;(9) 
CV 


A,u(9) 
Aj10 
d;10 
K\10 
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259.916 
] 
28 
9 


659.70 
1635.30 
0.29 
0.009 
0.033 


258.706 
] 
28 
6 


665.70 
1657.30 
0.29 
0.0094 
0.033 


297-912 
/ 
28 
4 


669.70 
1681.30 
0.28 
0.009 
0.032 


256.372 
] 
28 
8 


278.927 
6 
30 
g 


472.2] 
803.40 
0.37 
0.014 
0.037 


277.659 
6 
30 
10 


482.27 
823.40 
0.37 
0.0134 
0.036 


277.023 
6 
30 
1 


483.27 
852.40 
0.36 
0.013 
0.036 


271.363 
10 
50 

8 


136.481 
2 
12 
2 


135.94 
415.81 
0.25 
0.018 
0.074 


135.506 
2 
12 
1 


136.94 
426.81 
0.24 
0.0180 
0.074 


133.598 
2 
12 
3 


139.94 
435.81 
0.24 
0.018 
0.073 


133.679 
Z 
12 
1 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Eleventh Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Twelfth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


c,(10) 
Bj(10) 
1(10) 
5;(10) 
CV 


Aju, 10) 


A,11 
di11 
K;11 


a(11) 
Bd) 
4,11) 
oj(11) 
CV 


Aww 1) 


A\12 
12 
K;12 


a,(12) 
Bj)(2) 
H(12) 
9,(12) 
Cy 
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677.70 
1701.30 
0.28 
0.009 
0.032 


256.381 
7 
28 
5 


682.70 
1724.30 
0.28 
0.009 
0.032 


255.268 
7, 
28 
6 


688.70 
1746.30 
0.28 
0.00913 
0.032 


491.27 
894.40 
0.35 
0.013 
0.036 


265.901 
10 
50 

8 


499.27 
936.40 
0.35 
0.013 
0.036 


260.82 


509.27 
976.40 
0.34 
0.01231 
0.036 


140.94 
446.81 
0.24 
0.018 
0.073 


131.885 
2 
12 
1 


141.94 
457.81 
0.24 
0.017 
0.101 


142.94 
468.81 
0.23 
0.01709 
0.101 


OBJECTIVE: MAXIMIZE NUMBER OF ILLEGAL SHIPS FOUND 


Total Assets (n) 
Real X; (j=1,2,3) 
Variable description 
Area information 


Ship arrival rate 
Assets inspection capacity 


Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 


Initial Beta parameters 

Alpha 

Beta 

First Inspection 

Illegal ship rate 

Initial assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 

Updated illegal rate 

Second Inspection 

New assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


0.2 


95% ODDS RATIO BOUNDS 
LOWER 
85 


Cap, 
14,(0) 
6;(0) 


a;(0) 
B,(9) 


A,u,(0) 
A(1) 

dj(1) 

E(Kj(1)] 

Kj(1) 
P(X>p,(0))/P(XJ<p,(0)) 
H(1) 

o,(1) 


aj(1) 
Bi) 
Ajual 1 ) 


A2 
42 
E{Kj(2)} 

K2 
P(X>4(1))/PO9<4 (1) 
H;(2) 

0,(2) 


a,(2) 
B;(2) 
ju,(2) 


Area 1 
0.42 
2.92 


Area 1 


300 


4 
0.28 
0.028 


71.72 
184.42 


252 


28 
3.6 


0.81 
0.28 
0.03 


71.72 
184.42 
252 


28 
5.6 


1.24 
0.28 
0.028 


71.72 
184.42 
252 
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Area 2 
0.17 
1.17 


Area 2 


750 


o 
0.34 
0.034 


65.66 
127.46 


255 


30 
6.0 


0.47 
0.34 
0.034 


65.66 
127.46 
255 


30 
6.0 


0.59 
0.34 
0.034 


65.66 
127.46 
255 


Area 3 
0.83 
3.52 


Area 3 


550 


6 
0.23 
0.023 


76.77 
257.01 


127 


12 
2.4 


az 
0.23 
0.023 


76.77 
257.01 
127 


12 
2.4 


152 
0.23 
0.019 


76.77 
257.01 
127 


Third Inspection 


New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


Fourth Inspection 


New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


Fifth Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


A3 
d3 
ELKi(3) 

K(3) 
P(Xj>H,(2))/P(X9<H,(2)) 
(3) 

G;(3) 


a(3) 
B;(3) 


Aj4) 
dic4) 
E[Kj(4)} 

Ki4) 
P(Xj>n,(3))/P(X)<1,(3)) 
(4) 

G,(4) 


a(4) 
B;(4) 
Anil) 


Ai) 
dj9) 
E(Kj(5)] 

Ki) 
P(X>1,(4))/P(X)<p,(4)) 
H;(5) 

G,(5) 


a;(5) 
B;(5) 
Ajpi{9} 


28 
5.6 


0.81 
0.28 
0.028 


71.72 
184.42 
252 


28 
5.6 


0.53 


0.28 
0.03 


We72 
184.42 
252 


5 
20 
4.0 

4 

0.69 


0.28 
0.03 


71.72 
184.42 
252 
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30 
6.0 


0.47 
0.34 
0.034 


65.66 
127.46 
255 


30 
6.0 


0.47 


0.34 
0.03 


65.66 
127.46 
255 


20 
40 


0.31 
0.34 
0.03 


65.66 
127.46 
255 


12 
2.4 


0.83 
0.23 
0.023 


76.77 
257.01 
127 


12 

2.4 

10 
4.39 
0.83 
0.42 


0.17 
-0.03 
458 


36 
14.4 


0.00 
0.36 
0.18 


2.19 
3.88 
199 


—_ ganh 


mo on OD oO &S WW PO 


ami eam 
—_ C9 


—" 
Ro 


Do a nay Oot WO Ht | OS 


ob eh 
no =— & 


CV Area 1 
0.033 
0.033 
0.039 
0.033 
0.033 
0.033 
0.033 
0.033 


0.032 
0.032 


0.032 
0.032 


Mean Area 1 
0.3 
0.30 
0.3 
0.29 


0.29 
0.29 


0.29 
0.29 
0.29 
0.28 
0.28 
0.28 
0.28 


CV REDUCTION 


CV Area 2 


0.037 
0.037 
0.046 
0.037 
0.037 
0.037 
0.037 
0.036 


0.036 
0.036 


0.036 
0.036 


MEAN 


Mean Area 2 


0.4 
0.40 
0.39 
0.39 


0.38 
0.38 


0.37 
0.37 
0.37 
0.36 
0.35 
0.35 
0.28 


CV Area 3 
0.079 
0.078 
0.088 
0.076 
0.075 
0.075 
0.074 
0.074 


0.073 
0.073 


0.101 
0.101 


Mean Area 3 
0.25 
0.25 
0.25 
0.25 


0.25 
0.25 


0.25 
0.25 
0.24 
0.24 
0.24 
0.24 
0.23 
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INSPECTIONS MISSIONS RESULTS SCENARIO: 


ASSETS 
INSPECTION Areal Area 2 Area 3 
IST 7.0 6.0 2.0 
2ND 70 6.0 2.0 
3RD 7.0 6.0 2.0 
4TH 7.0 6.0 2.0 
STH 5.0 4.0 6.0 
POSTERIOR ODDS RATIO 
INSPECTION Area 1 Area 2 Area3 
IST 0.812 0.469 1.523 
2ND 1.239 0.592 1.523 
3RD 0.81 0.47 0.83 
4TH 0.529 0.469 4.389 
STH 0.69 0.31 0.00 
MEAN 
INSPECTION Area 1 Area 2 Area 3 
toh 0.280 0.340 0.230 
2ND 0.280 0.340 0.230 
3RD 0.280 0.340 0.230 
4TH 0.280 0.340 0.833 
STH 0.280 0.340 0.361 


EXPECTED NUMBER OF ILLEGAL SHIPS 


INSPECTION Area I Area 2 Area 3 Total 
IST 7 7 5 19 
2ND 9 8 5 22 
3RD 6 6 2 19 
4TH 6 6 2 14 
STH 4 4 14 22 

96 
# ILLEGAL SHIPS FOUND 

INSPECTION Area 1 Area 2 Area 3 Total 
IST 7 7 5 19 
2ND 9 8 5 22 
3RD 7 7 2 16 
4TH 5 7 10 22 
STH 4 2 13 19 

98 
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OBJECTIVE: REDUCE UNCERTAINTY OF X ; 


Total Assets(n) 20 (p Value 0.1 
Real X; Q=1,2,3) 0.2 

Variable description Symbol Areal Area2 Area3 
Area information 

Ship arrival rate i 900 = 750 550 

Assets inspection capacity Cap, 4 5 6 

Initial mean fraction of illegal ships 14,(0) 0.3 0.4 0.25 

Initial S.D. fraction of illegal ships 6;(0) 0.010 0.015 0.020 
Initial Beta parameters 

Alpha a;(0) 629.70 426.27 116.94 

Beta Bi(0) 1469.30 639.40 350.81 
First Inspection 

Illegal ship rate AjH,(0) 270 300 137.5 

Initial assets allocation A; 9 8 3 

# of ships to be inspected dj1 36 40 18 

# of illegal ships found i 3 8 8 
Updated Beta distribution 

Alpha (1) 632.70 434.27 124.94 

Beta Bi) 1502.30 671.40 360.81 

Updated mean fraction of illegal ships (1) 0.30 0.39 0.26 

Updated S.D. fraction of illegal ships 6;(1) 0.010 0.015 0.020 

Coefficient of variation Cy 0.033 0.037 0.077 

Updated illegal rate AjH,(1) 266.712 295 141 
Second Inspection 

Assets allocation Aj2 9 8 3 

# of ships to be inspectcd dj2 36 40 18 

# of illegal ships found Kj2 7 5 3 
Updated Beta Distribution 

Alpha o;(2) 639.70 439.27 127.94 

Beta Bi(2) 1531.30 706.40 375.81 

Updated mean fraction of illegal ships U;(2) 0.29 038 0.25 

Updated S.D. of fraction of illegal ships 6;(2) 0.010 0.014 0.019 

Coefficient of variation CV 0.033 0.037 0.076 
Third Inspection 

New illegal rate Ap(2) 261 285 138 

Assets allocation Aj3 g 8 3 

# of ships to bec inspected dj3 36 40 18 

# of illegal ships found Kj(3) 6 13 1 
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Updated Beta distribution 


Alpha a;(3) 645.70 452.27 128.94 
Beta B,(3) 1561.30 733.40 392.81 
Updated mean fraction of illegal ships L;(3) 0.29 0.38 0.25 

Updated S.D.fraction of illegal ships G,(3) 0.010 0014 0.019 


Coefficient of variation CV 0.039 0.046 0.086 
Fourth Inspection 
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New illegal rate Au(3) 261 285 =: 137.5 

Assets allocation Aj4 9 8 3 

# of ships to be inspected di4 36 40 18 

# of illegal ships found Kj4 3 6 4 
Updated Beta Distribution 

Alpha ot;(4) 648.70 458.27 132.94 

Beta B;(4) 1594.30 767.40 406.81 

Updated mean fraction of illegal ships 14;(4) 0.29 0.37 0.25 

Updated S.D. fraction of illegal ships G;(4) 0.010 0014 0.019 

Coefficient of variation CV 0.033 0.037 0.075 
Fifth Inspection 

New illegal rate Au(4) 261 77a) | (AGT 

Assets allocation Ajo 9 8 3 

# of ships to be inspected djs 36 40 18 

# of illegal ships found Kj9 10 4 ] 
Updated Beta Distribution 

Alpha o;(5) 658.70 462.27 139.94 

Beta Bj(S) 1620.30 803.40 417.81 

Updated mean fraction of illegal ships u(5) 0.29 0.37 0.25 

Updated S.D. fraction of illegal ships G,(5) 0.0095 0.0135 0.0183 

Coefficient of variation CV 0.033 0.037 0.073 
Sixth Inspection 

New illegal rate A(5)  —- 260.127 273.927 137.993 

Assets allocation A6 9 8 3 

# of ships to be inspected dj6 36 40 18 

# of illegal ships found K;6 ] 10 4 
Updated Beta Distribution 

Alpha a;(6) 665.70 472.27 143.94 

Beta Bi(6) 1649.30 833.40 431.81 

Updated mean fraction of illegal ships 14,(6) 0.29 0.36 0.25 

Updated S.D. fraction of illegal ships G;(6) 0.009 0013 0.018 

Coefficient of variation CV 0.033 0.037 0.072 


Seventh Inspection 
New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Eigth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Ninth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 
# of illegal ships found 
Updated Beta Distribution 
Alpha 
Beta 
Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Coefficient of variation 
Tenth Inspection 


New illegal rate 
Assets allocation 
# of ships to be inspected 


# of illegal ships found 
Updated Beta Distribution 


Aip;(6) 
Aj/ 
dj7 
Kj/ 


a;(7) 
Bi(7) 
H(7) 
5,(7) 
Cy 


Aipi{7) 
A8 
d.8 
K8 


a;(8) 
B,(8) 
H;(8) 
9,(8) 
oN/ 


Aiu,(8) 
Aj9 
dd 
K9 


a;(9) 
Bi(9) 
H,(9) 
G;(9) 
ey 


Aiu(9) 
Aj10 
d;10 
K;10 
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258.803 
9 
36 
6 


671.70 
1679.30 
0.29 
0.009 
0.033 


257.137 
9 
36 
6 


677.70 
1709.30 
0.28 
0.0092 
0.032 


255.522 
9 
36 
8 


685.70 
1737.30 
0.28 
0.009 
0.032 


254.697 
9 
36 
] 


271.279 
8 
40 
6 


478.27 
867.40 
0.36 
0.013 
0.037 


266.559 
8 
40 
8 


486.27 
899.40 
0.35 
0.0128 
0.037 


263.195 
8 
40 
7 


493.27 
932.40 
0.35 
0.013 
0.036 


259.493 
10 
50 

8 


146.94 
446.81 
0.25 
0.018 
0.072 


136.111 
3 
18 
5 


151.94 
459.81 
0.25 
0.0175 
0.070 


136.601 
3 
18 
6 


157.94 
471.81 
0.25 
0.017 
0.068 


137.937 
3 
18 
4 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Eleventh Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


Twelfth Inspection 


New illegal rate 

Assets allocation 

# of ships to be inspected 
# of illegal ships found 


Updated Beta Distribution 


Alpha 

Beta 

Updated fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Coefficient of variation 


a(1 0) 
B;(10) 
H,(10) 
o;(10) 
CV 


2i11,(10) 


As11 
11 
K11 


ay(11) 
Bi1) 
(11) 
o,(11) 
CV 


Aw(11) 
Aj12 
d12 
K12 


a(12) 
B12) 
H(12) 
o,(12) 
CV 
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692.70 
1766.30 
0.28 
0.009 
0.032 


253.53 
9 
36 
iB 


703.70 
1791.30 
0.28 
0.009 
0.032 


253.84 


36 


712.70 
1818.30 
0.28 
0.00894 
0.032 


901.27 
974.40 
0.34 
0.012 
0.036 


254.766 
10 
90 

8 


909.27 
1016.40 
0.33 
0.012 
0.036 


250.35 
10 
50 
14 


923.27 
1052.40 
0.33 
0.01186 
0.036 


161.94 
485.81 
0.25 
0.017 
0.068 


137.5 


18 


165.94 
499.81 
0.25 
0.017 
0.101 


137.087 
3 
18 
4 


169.94 
913.81 
0.25 
0.01652 
0.101 


OBJECTIVE: MAXIMIZE NUMBER OF ILLEGAL SHIPS FOUND 


Total Assets (n) 
Real Xj §=1,2,3) 
Variable description 
Area information 


Ship arrival rate 
Assets inspection capacity 


Initial mean fraction of illegal ships 
Initial S.D. fraction of illegal ships 


Initial Beta parameters 

Alpha 

Beta 

First Inspeetion 

Illegal ship rate 

Initial assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 

Updated illegal rate 

Seeond Inspeetion 

New assets allocation 

# of ships to be inspected 

Expected number of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


95% ODDS RATIO BOUNDS 


LOWER 
95 


Symbol 


dj 
Cap, 
H1;(0) 
3;(0) 


a(0) 
B;(0) 


AjH,(0) 
A(1) 
di(1) 
E(Kj(1)] 
K;(1) 


P(X}>H;(0))/P(XJ<H,(0)) 


H(1) 
o,(1) 


aj(1) 
By.) 


Aj2 
2 
E[Kj(2)] 
K2 


P(Xj>45(1))/P(CXJ<H(1)) 


H(2) 
or 2) 


a,(2) 
By(2) 


Area | 
0.26 
1.78 


Area | 


900 


4 
0.28 
0.028 


71.72 
184.42 


252 


36 
7.2 


0.51 
0.28 
0.03 


yes 2 
184.42 
252 


36 
ye 


0.77 
0.28 
0.028 


71.72 
184.42 
252 
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Area 2 
0.07 
0.55 


Area 2 


750 


5 
0.34 
0.034 


65.66 
127.46 


255 


40 
8.0 


0.21 
0.34 
0.034 


65.66 
127.46 
255 


40 
8.0 


0.27 
0.34 
0.034 


65.66 
127.46 
200 


Area 3 
0.62 
2.60 


Area 3 


550 


6 
0.23 
0.023 


76.77 
257.01 


127 


3.6 


1.15 
0.23 
0.023 


76.77 
257.01 
127 


3.6 


Tan 
0.23 
0.019 


76.77 
257.01 
127 


Third Inspection 


New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 
Updated Beta parameters 

Alpha 

Beta 


New illegal rate 


Fourth Inspection 


New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


Fifth Inspection 

New assets allocation 

# of ships to be inspected 

Expected value of illegal ships 

# of illegal ships found 

Posterior Odds Ratio 

Updated mean fraction of illegal ships 
Updated S.D. fraction of illegal ships 


Update Beta Distribution 
Alpha 
Beta 


New illegal rate 


A3 
d3 
E[Kj(3)] 

K;(3) 
P(X}>n;(2))/P(X}<H,(2)) 
H;(3) 

6;(3) 


a(3) 
B;(3) 
Au(3) 


Aj) 
dic4) 
E[Kj(4)] 

Ki4) 
P(X>14,(3))/P(XJ<H,3)) 
H;(4) 

0;(4) 


a;(4) 
B;(4) 
Au) 


Aj) 
dio) 
E[Kj(5)} 

KS) 
P(X}>4,(4))/P(X3<H,(4)) 
H(5) 

6,(5) 


a;(5) 
B;(5) 
Aiu,(5) 


36 
daa 


0.51 
0.28 
0.028 


71.72 
184.42 
252 


36 
7.2 


0.33 


0.28 
0.03 


VIZ 
184.42 
252 


32 
6.4 


0.52 
0.28 
0.03 


71.72 
184.42 
252 
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40 
8.0 


0.21 
0.34 
0.034 


65.66 
127.46 
255 


40 
8.0 


0.21 


0.34 
0.03 


65.66 
127.46 
255 


30 
6.0 


0.23 
0.34 
0.03 


65.66 
127.46 
255 


3.6 


0.62 
0.23 
0.023 


76.77 
257.01 
127 


3.6 
10 
3.21 
0.56 
0.28 


1.22 
0.98 


306 


36 
14.4 


0.01 
0.36 
0.18 


2.19 
3.88 
199 


—_— gmp 


ao ao un OO & WF PMY 


a_i ooh 
—_— CC 


a" 
INO 


Co Oo NN OO Oh WH NH = Oo 


— kek 
LS | 


CV Area 1 
0.033 
0.033 
0.039 
0.033 
0.033 
0.033 
0.033 
0.032 


0.032 
0.032 


0.032 
0.032 


Mean Area 1 
0.3 
0.30 
0.29 
0.29 


0.29 
0.29 


0.29 
0.29 
0.28 
0.28 
0.28 
0.28 
0.28 


CV REDUCTION 


CV Area 2 


0.037 
0.037 
0.046 
0.037 
0.037 
0.037 
0.037 
0.037 


0.036 
0.036 


0.036 
0.036 


MEAN 


Mean Area 2 


0.4 
0.39 
0.38 
0.38 


0.37 
0.37 


0.36 
0.36 
0.35 
0.35 
0.34 
0.33 
0.28 


CV Area 3 
0.077 
0.076 
0.086 
0.075 
0.073 
0.072 
0.072 
0.070 


0.068 
0.068 


0.101 
0.101 


Mean Area 3 
0.25 
0.26 
0.25 
0.25 


0.25 
0.25 


0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
0.25 
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INSPECTIONS MISSIONS RESULTS SCENARIO: 


ASSETS 
INSPECTION Area 1 Area 2 Area 3 
1ST 9.0 8.0 3.0 
2ND 9.0 8.0 3.0 
3RD 9.0 8.0 3.0 
4TH 9.0 8.0 3.0 
STH 8.0 6.0 6.0 
POSTERIOR ODDS RATIO 
INSPECTION Area 1 Area 2 Area3 
1ST 0.507 0.212 1.148 
2ND 0.774 0.271 1.148 
3RD 0.51 0.21 0.62 
4TH 0.327 0.212 3.213 
STH 0.52 0.23 0.01 
MEAN 
INSPECTION Areal Area 2 Area 3 
ISeé 0.280 0.340 0.230 
2ND 0.280 0.340 0.230 
3RD 0.280 0.340 0.230 
4TH 0.280 0.340 0.556 
STH 0.280 0.340 0.361 


EXPECTED NUMBER OF ILLEGAL SHIPS 


INSPECTION Area I Area 2 Area 3 Total 
IST 7 7 5 19 
2ND g 8 5 22 
3RD 7 8 4 19 
4TH ii 8 4 19 
STH 6 6 14 27 

106 
# ILLEGAL SHIPS FOUND 

INSPECTION ‘Areal Area2 ‘Area3 Total 
IST 7 7 5 19 
2ND 9 8 5 22 
3RD 7 7 2 16 
4TH 5 7 10 22 
STH 6 4 13 23 

102 
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APPENDIX B. CUMULATIVE RIGHT TAIL BINOMIAL PROBABILITIES 


dj(t+1) = 12 








| |_| Cumulative Right Tail Binomial Probabilities | = = 

en ee 
ewes 

_ eae eee eee eee 
PT ft tates | 
Lo ae eee eee eee 

| 1 | 02s | 0023 | 0002 | 0.000 | 0.000 | 000 | 0. 
| 1 | 037 | 0.081 | oor | 001 | 0.000 || 0.000 | 0.000 | 
| tT 0524 | go | 0032 | 000s | 000 | 0.000 *|| 0.000 
| 1 | 0632 | 0249 | oes | iz | 02 | 000 | 0.000 
P| oe | 0341 | ont | 0.026 | 0004 | .001 | 0.000 
| tf oes | 04s) | oe? | 0.048 |S o9_— | 0.001 | 0.000_— 
| iT eae | siz | 230 | Sos | ore | 0.003, |= 0.000_— 
| | 0877 | ses | zee | oti) =| 0.031 | 0.006 =| ool 
| 1 | ogee | eee | sro | ss | 0.048 | 02 | 0.002 
| i | oti | 0726 | 0az | 20s | 0073 | 001d | 0m 
| it | oesg | ove | ost | 261 | o.t02 | 0.030 | 0.007 
| 1 | ages | 0822 | ose | 0320 | se | 0s | oN 
| 1 | 0973 | 088s | eso | ose2_ | 079 | 008s | 08 
1 {| 0981 | 0890 | 0898 | 0445 | 0.225 | 0.089 
tf oess fogs | 747 | som | 0276 | 08 
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